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FOREWORD 


This  is  the  final  report  covering  the  technical  work 
performed  by  the  Department  of  Aerospace  Engineering  and 
Applied  Mechanics  under  Contract  No.  DAAG-69-C-0016  (U.S. 
Research  Office  Durham) . The  contract  was  initiated  under 
Project  Themis  No.  DAHC04-69-C-0016 , entitled  "Internal 
Aerodynamics,  Heat  Transfer  and  High  Temperature  Materials 
in  Air-Breathing  Engines.”  The  work  was  supervised  by 
Dr.  Widen  Tabakoff,  Professor  of  the  Department  of  Aerospace 
Engineering  and  Applied  Mechanics.  The  report  covers 
work  conducted  from  September  20,  1968  to  June  30,  1976. 

The  contents  of  the  technical  reports  written  during 
the  contract  period  have  been  incorporated  in  this  report 
through  the  provision  of  abstracts.  In  addition,  the 
technical  publications  in  scientific  journals  are  given 
with  the  first  page  of  the  corresponding  publication. 

The  scientific  and  technical  personnel  participating 
in  this  research  project  are  listed  in  Appendix  A,  and 
the  students  obtaining  degrees  in  Appendix  B. 
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A brief  review  of  the  Project  Themis  Research  Program  in 
Air  Breathing  Propulsion  at  the  University  of  Cincinnati  is 
presented.  The  review  covers  the  period  from  September  1968 
to  June  1976.  The  research  presented  herein  is  concerned 
with  improving  the  performance  of  air-breathing  propulsion 
systems.  Theoretical  and  experimental  studies  of  basic 
phenomena  upon  which  the  performance  is  dependent,  were 
carried  out.  These  studies  covered  the  following  areas: 

1.  Particulate  Flows  in  Propulsion  Systems 

2.  Cascade  Aerodynamics 

3.  Jet  Mixing  Flow 

4.  High  Temperature  Materials 

5.  Gas  Turbine  Blades  and  Combustion  Chambers 
Heat  Transfer. 
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PARTICULATE  FLOWS  IN  PROPULSION  SYSTEMS 


I . 

The  study  of  gas  particle  flew  in  turbomachines  is  of 
great  interest.  As  a result  of  the  development  of  high  energy 
propellants  for  propulsion  systems,  the  products  of  combustion 
may  partially  consist  of  finely  divided  particles.  For 
example,  solid  propellant  rockets  are  sometimes  used  as  starting 
devices  for  propulsion  systems  employing  turbines.  Air- 
breathing  engines  operating  in  desert  areas,  where  the  inlet 
air  flow  may  contain  sand  particles,  provide  another  example. 

Most  of  these  types  of  flows  involve  changes  in  the  gas 
velocity  and  temperature.  Gas  particle  interaction,  through 
viscous  drag  and  heat  transfer  produces  corresponding  changes 
in  the  properties  of  particles.  These  particles  are  swept 
through  the  turbine  by  the  gas  flow,  lagging  behind  the  gas 
in  velocity  and  temperature.  Both  these  effects  cause  a 
deterioration  in  the  performance  of  turbomachines.  Furthermore, 
these  particles  may  cause  blade  erosion,  and  consequently, 
failure.  Even  if  the  erosion  is  not  severe  enough  to  cause 
rupture,  changes  in  blade  geometry  are  of  concern,  due  to  their 
influence  on  performance. 

The  study  of  the  trajectories  and  velocities  of  solid 
particles  suspended  in  a fluid  flow  through  cascades  and  blade 
rows  is  necessary  to  investigate  the  erosion  damage  sustained 
by  the  blades.  In  general,  the  trajectories  and  velocities  of 
the  particles  depend  upon  the  shape  of  the  cascade,  the  particle 
and  flow  inlet  conditions,  the  flow  angle  of  attack,  the 
particle  material  density,  the  mean  diameter  and  the  location  of 
the  initial  particle  collision.  The  particles  are  more  likely 
to  follow  the  fluid  streamlines  when  their  material  density  is 
of  the  same  order  of  magnitude  as  that  of  the  fluid  and  when 
their  mean  diameter  is  small.  The  pressure  distribution  on 
the  blade  surfaces  changes  when  solid  particles  are  introduced 
into  the  gas  stream.  Consequently,  the  performance  of  the  turbine 
or  compressor  in  an  air-breathing  engine  will  change  as  well. 
Knowing  the  effect  of  the  particles  on  the  blade  pressure 
distribution  will  be  very  useful  for  the  design  purposes. 

The  research  work  performed  in  the  above  mentioned  area 
is  reported  in  the  following  abstracts  and  references. 
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Hamed , A.  and  Tabakoff,  W.  , "Momentum  Integral  Theory  for  Gas 
Particle  Flow,"  Project  Themis  Report  No.  69-3,  September 
1969.  (AD  6 9 6142) 

Abstract 

A theoretical  treatment  of  the  boundary  layer  for  the  two- 
dimensional  gas  particle  flov;  is  presented.  The  momentum 
integral  equations  of  the  boundary  layer  with  incompressible 
gas  phase  are  derived,  and  a numerical  calculation  is  carried 
out  for  the  case  of  a flat  plate. 


Hussein,  M.  and  Tabakoff,  W. , " Th e Properties  of  a Gas  Particle 
Flow  in  Cascade,"  Project  Themis  Report  69-4,  September 
1969.  (AD  697164) 

Abstract 

An  approximate  method  is  developed  for  calculating  the 
flov;  properties  of  gas-particle  mixture  flowing  over  turbine 
blades  in  a cascade.  Momentum  and  heat  transfer  between  the 
gas  and  particles,  as  well  as  compressibility  effects  were 
considered.  An  investigation  was  made  of  the  effect  of 
particle  concentration,  mean  diameter,  particle  material 
density,  particle  inlet  velocity  and  temperature  on  the  gas 
particle  flow  properties  on  a turbine  blade. 


Tabakof f,  W.  and  Hussein,  M. , "An  Experimental  Study  of  the 
Effect  of  Solid  Particles  on  the  Pressure  at  Blade  Surfaces  in 
Cascade,"  Project  Themis  Report  No.  70-8,  March  1970.  (AD  703896) 

Abstract 

Data  representing  the  effect  of  solid  particle  flow  on  the 
pressure  distribution  over  blade  surfaces  in  cascade  are  presented. 
Various  inlet  mass  flow  concentrations  were  tested.  Different 
particle  sizes  were  used,  namely  50,  300  and  1000  microns  in 
order  to  determine  their  effect  on  the  pressure  distribution. 


Tabakoff,  W.  and  Hussein,  M.,  "Experimental  Investigation  of  the 
Trajectories  and  Velocities  of  Solid  Particles  Entrained  by 
Fluid  Flows  in  Cascade  Nozzles,"  Project  Themis  Report  No.  70-12, 
August  1970.  (AD  711121) 

Abstract 

An  experimental  investigation  was  made  to  determine  the 
trajectories  and  velocities  of  solid  particles  suspended  in  a 
fluid  flowing  through  a cascade  nozzle. 
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Hamed,  A.  and  Tabakoff,  V/.  , " The  B oundary  Layer  of  Particulate 

F lov:  in  Cascade  Nozzles,"  Project  Themis  Report  No.  70-16 , 

November  1970.  (ad  715721) 

Abstract 

The  boundary  layer  of  gas-particle  flows  in  field  of  changing 
pressures  is  analyzed.  The  analysis,  which  utilizes  the  momentum 
integral  method,  can  be  used  to  study  the  boundary  layer  of  any 
particulate  flow  with  incompressible  gas-phase  provided  the  flow 
properties  outside  the  boundary  layer  are  known.  Numerical 
integration  of  the  two  momentum  integral  equations  was  carried 
out  for  the  case  of  particulate  flow  in  cascades.  For  this 
example,  the  results  indicate  that  introduction  of  particles 
leads  to  an  increase  in  the  gas  boundary  layer  thickness  and  that 
the  gas  boundary  layer  characteristics  are  usually  more  sensitive 
to  particle  concentration  than  any  other  particulate  flow  para- 
meter. However,  the  particle  boundary  layer  characteristics  are 
more  sensitive  to  mean  particle  diameter  than  to  particle  con- 
centration. Other  parameters  such  as  inlet  particle  gas  velocity 
ratio  and  particle  solid  material  density  have  more  pronounced 
effects  on  the  particle  boundary  than  on  that  of  the  gas. 


Hussein , M.  and  Tabakoff,  W.  , "Gas-Particle  Suspension  Properties 
on  a Blade-to-Blade  Stream  Surface  of  a Cascade  Nozzle,”  Project 
Themis  Report  No.  70-17,  December  1970.  (AD  715971) 

Abstract 

The  gas-particle  flov;  properties,  namely  flov;  pressure,  gas 
and  particles  velocities  and  temperatures  were  calculated  theo- 
retically in  a blade  to  blade  stream  surface  of  a cascade  nozzle. 
The  effects  of  the  particle  concentration,  mean  diameter  and 
material  density  on  the  flow  properties  at  any  point  in  the 
cascade  nozzle  were  investigated. 


Tabakoff,  W. , Hamed,  A.  and  Hussein,  M.F.,  "Experimental  Investi- 
gation of  Gas  Particle  Flov;  in  a Simulated  Axial  Flow  Compressor 
Stage,"  Project  Themis  Report  No,  71-21,  July  1971^  (AD  728489) 

Abstract 

An  experimental  study  of  the  effect  of  particle  size  and  fluid 
velocity  on  the  trajectories  and  velocities  of  solid  particles 
entrained  by  a fluid,  and  passing  through  a simulated  axial  flow 
compressor  stage,  is  carried  out.  Conclusions  concerning  areas  of 
blade  surface  subjected  to  erosion  damage  due  to  particle  collisions 
are  drawn.  Pressure  measurements  of  the  particulate  gas  flow  on 
the  blade  surfaces  were  obtained  for  various  values  of  particle  con- 
centration for  two  different  positions  of  the  second  cascade  row. 

A subsonic  cascade  wind  tunnel  was  used  for  testing  the  particulate 
flov;  through  the  simulated  compressor  stage. 


PROPERTIES  AND  PARTICLE  TRAJECTORI 
FLOWS  IN  CASCADES 


OF  GAS-PARTICLE 


An  approximate  method  for  calculating 
the  flow  properties  of  gas-particle  mix- 
ture flowing  over  blades  in  a cascade  is 
developed.  The  momentum  and  heat  transfer 
between  the  gas  and  particles  as  well  as 
compressibility  effects  are  taken  into 
consideration.  The  effect  of  particle 
concentration,  mean  diameter,  particle 
material  density,  particle  inlet  velocity 
and  temperature  on  the  gas  paiticle  flow 
properties  are  investigated.  The  flow 
pressure  distribution  on  the  blade  surface 
and  the  particle  velocities  and  trajecto- 
ries in  the  nozzle  are  determined  experi- 
mentally. The  experiments  are  performed 
in  a wind  tunnel  and  in  a water  table. 

The  solid  particle  sizes  range  from  5 to 
1C00  microns  and  the  particle  concentra- 
tion range  from  0.03  to  0.2. 


W.  TAPAKOFF  and  M.F.  HUSSEIN 
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Trajectories  oi’  Particles  Suspended  in 
Fluid  Finn  Through  Cascades 

W.  Tahakoi  k*  and  M.  F.  HussurN't 
University  oj  Cincinnati,  Cincinnati,  Ohio 

Introduction 

HU  HE  study  oi  the  trajectories  and  velocities  of  solid  par- 
tides  suspended  in  a fin  id  flow  passing  through  a cascade 
nozzle  is  of  importance  to  the  investigation  of  erosion  damage 
sustained  by  the  blades.  In  general,  the  trajectories  and 
velocities  of  the  particles  depend  upon  the  slope  of  the  cas- 
cade nozzle,  particle  and  flow  ir.iet  conditions,  particle  ma- 
teria! density,  mean  diameter,  angle  of  attack  and  initial 
place  of  collision.  The  particles  are  more  likely  to  follow  the 
fluid  streamlines  when  their  material  density  is  of  the  same 
order  of  magnitude  as  that  of  the  fluid  and  when  their  mean 
diameter  is  small. 

In  this  experimental  investigation,  a cascade  row  of  turbine 
blades  was  mounted  in  a cascade  tunnel  to  produce  the  desired 
gas-particle  flow.  A high-speed  camera  was  used  to  photo- 
graph the  flow.  The  film  analysis  provided  the  data  for  the 
particle  velocities  and  the  particle  paths  through  the  cascade 
nozzle. 

Test  Facilities  and  Conditions 

A special  subsonic  cascade  wind  tunnel,  designed  for  gas- 
particle  studies,  was  used  for  this  scries  of  tests.  Reference 
1 gives  a detailed  description  of  the  tunnel  and  the  airfoil 
used  in  the  cascade  nozzle. 

A high  speed  motion  picture  camera  with  a maximum  speed 
of  6000  frames-sec  at  2'20-v  input  was  used  to  photo- 
graph the  particulate  flow.  Solid  particles  of  1000-fi  mean 


diameter  ai  d 1.1-g/cm*  density  were  used.  The  range  of 
particle  diameters  available  was  limited  by  the  Ret  that  if 
particles  of  smaller  diameters  were  used,  it  was  impossible 
to  determine  their  trajectories  from  a film  analysis.  A direct 
lighting  technique  was  used  in  the  experiments,  and  the 


Fig.  1 Trajectories  of  particles  entrained  by  gas  flow 
(patli  a). 


Presented  as  Paper  70-712  at  the  AIAA  0th  Propulsion 
Joint  Specialist  Conference,  San  Diego,  Calif.,  June  15-10,  1070; 
submitted  July  13,  1970.  This  work  was  sponsored  under  Proj- 
ect Themis  Contract  Number  DAHCO  4-69C-0010,  U.S.  Array 
Research  Office — Durham. 
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Measurement:;  of  P.irticulatecl  Gus 
Flow  Pressure  on  Cascade  Nozzles 

W.  Taiiakoit*  and  M.  F.  Hussein] 
University  0/  Cincinnati,  Cincinnati,  Ohio 

Cascade  I unucl  and  Model  Description 

A SPECIAL  subsonic  cascade  tunnel  was  built  which  in- 
— - corporates  a device  for  injection  of  solid  particles  (Fig. 
1).  The  cascade  dimension : and  the  pressure  probe  locations 
arc  shown  in  Fig.  2.  Further  test  facility  information  may  be 
obtained  from  Kef.  1. 

Instrumentation 

Manometer  readings  were  re  orded  by  camera  after  a 
steady  state  was  reached.  The  primary  mass  flow  was  mea- 
sured by  an  orifice  meter  located  ahead'  of  the  settling  cham- 
ber. The  primary  flow  temperature  was  measured  with  a 
standard  thermocouple  located  in  the  settling  chamber.  An 
electronic  counter  was  employed  to  record  the  time  of  par- 


Prosented  as  Paper  70-7i2  at  the  AT  A A f.th  Propulsion  Joint 
Specialist  Conference,  San  Diego,  Calif.,  Juno  15-19,  1970; 
submitted  July  13,  1970;  revision  received  October  13,  1970. 
This  work  was  sponsored  under  Project  Themis  Contract  No 
DAIJC-01-69C-001G,  U.S.  Army  Research  Office— Durham. 

* Piofessor,  Department  of  Aerospace  Engineering.  Associate 
Fellow  AIAA. 

t Graduate  Research  Assistant,  Department  of  Aerospace 
Engineering.  Student  Member  AIAA. 
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SIMULATION  OF  ENVIRONMENTAL  SOLID-PARTICLES 
TRAJECTORIES  AND  VELOCITIES  THROUGH  AN  AXIAL  FLOW 
COMPRESSOR  STAGE,  AND  THE  PRESSURE  DISTRIBUTION  ON  BLADES 

Abstract 

An  experimental  investigation  of  the  trajectories  and 
velocities  of  solid  particles  suspended  in  a fluid  passing 
through  an  axial  flow  compressor  stage  was  performed.  Such 
an  investigation  is  of  importance  to  the  study  of  erosion 
damage  sustained  by  the  blade.  Two  test  facilities  were  used 
for  this  sbudy:  a subsonic  cascade  wind  tunnel  for  compressible 

flow  and  a water  table  for  incompressible  flow.  From  the  test 
technique  it  may  be  concluded  that  the  present  existing 
theoretical  analysis  for  particle  trajectories  through  a 
compressor  stage  is  not  valid.  The  wind  tunnel  test  simulation 
is  much  better  than  the  water  table  and  may  be  used  for 
predicting  particle  trajectories. 


BY 

W,  TABAKOFF,  A,  HAMED  AND  M.F,  HUSSEIN 
UNIVERSITY  OF  CINCINNATI,  CINCINNATI,  OHIO 


MAY  18-20,  1971  TRENTON,  NEW  JERSEY 


Abstract 


A theoretical  method  was  developed 
for  predicting  the  pressure  distribution 
over  a blade  in  cascades  for  both  incom- 
pressible and  compressible  flow  with 
particles.  Experimental  results  were 
obtained  from  a cascade  wind  tunnel  equip- 
ped with  a solid  particle  injecting  system. 
Good  agreement  was  noted  between  the 
theoretical  and  experimental  pressure  dis- 
tribution for  both  incompressible  and  com- 
pressible flow  conditions.  This  agreement 
indicates  that  the  pressure  distribution 
over  a blade  surface  in  essence  nozzles 
increases  with  the  introduction  of  solid 
particles  in  the  flow.  The  increase  in 
pressure  due  to  the  particles  gives  rise 
to  forces  on  the  blades  indicating  that 
this  information  is  of  importance  to  the 
designer. 


THE  COMPARISON  BETWEEN  THEORETICAL  AND  EXPERIMENTAL 
PRESSURE  DISTRIBUTION  ON  A BLADE  IN  CASCADE  NOZZLE 
FOR  A PARTICULATE  GAS  FLOW 
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Effect  of  Suspended  Solid  Particles 
on  the  Properties  in  Cascade  Flow 

W.  TaHAKOFF*  and  M.  F.  1 1 UK.-JEfN  f 
University  of  Cincinnati , Cincinnati,  Ohio 

An  approximate  method  for  ealciilatin^  tlie  flow  properties  of  ga$-partiele  mixture  flowing 
over  hlades  in  a cascade  is  developed.  The  momentum  and  heat  t< ansfer  hetw ecu  the  gas  and 
particles  as  well  us  compressibility  effects  arc  taken  into  consideration.  The  effect  of  particle 
concentration,  mean  diameter,  particle  material  density,  particle  inlet  velocity  and  tcinpcia- 
t tire  on  the  gas  particle  flow  properties  are  inv«  s ligated. 


! 
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I 
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Nomenclature 

A ~ cross-sectional  area  of  stream  tube  (ft1) 

A'  — nonJimehskmul  cross-sectional  area  of  stream  tube 
o = the  ratio  of  tire  mass  flow  rate  of  jiartii Irs  to  the  total 
mass  (low  rate  of  the  gas  and  particle  mixture 
A,  - moan  particle  projected  area  (ft1) 

A.  mean  particle  surface  area  (ft8) 

Cd  - particles  drag  coeflicieut 

C„  - specific  heat  of  solid  particles  material  (I5tu/lb°li) 

Cve  — gas  specific  beat  at  con.  taut  pressure  (Htu/Ib°ll) 

Is  distance  along  airfoil  contour  between  two  successive 
points  (ft) 

dt  — time  in  which  the  particle  travels  a distance  5 (sec) 
dp  — particle  mean  diameter  (ft) 

F coefficient 

/ correction  factor  to  Stokes  drag  formula 
9 = gravity  constant 

h - coefficient  of  heat  transfer  between  the  particles  and  the 
gas  (Btu/lir  f ts!>Ii) 

J — mechanical  equivalent  of  heat 

k , ~ coefficient  of  conductivity  for  tfie  gas  (lltu/hr  ft°l!) 

m mass  of  one  particle  (lb) 

A'u  — Nusselt  number 

/’  «•  gas  particle  suspension  pressure  (lb/fl1) 

/"  = nondimensional  gas  particle  suspension  pressure 

p — pressure  of  gas-only  flow  (lb/fl1) 

p’  - nondimensional  gas-only  flow  pressure 

Pr  = Prandtl  number 

U j — gas  constant 

l!r  ---  Reynold.-,  number 

p = the  gas-only  flow  density  (lb/fl*) 

p'  = the  gas-only  flow  nondimensional  density 

P,  = the  gas  density  in  gas  particle  flow  (lb/ft*) 

Pp  — the  particles  density  (lb/fl1) 

Pp'  — the  nondimensional  particles  density 
Pp  = solid  particle  material  density  (lb/ft3) 
a «»  distance  along  airfoil  contour  (ft) 

J\  --  total  temperature  of  gas  (°R) 

T "■  the  gas-only  flow  tempo  attire  (°I!) 

T - the  nondimensional  temperature  of  gas-only  flow 
T , — the  gas  temperature  in  gas  particle  flow  (°K ) 

T./  = the  nondimensional  gas  temperature  in  gas  particle  flow 
Tp  = particle  temperature  (°I1) 
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Tp'  3=  nondimensional  particle  temperature 

t 3=  time  (see) 

u — gas-only  flow  velocity,  fps 

u'  — gus-only  flow  nondimensional  velocity 

ua'  — nondimensional  gas  velocity  in  gas  particle  flow 

v,  - gas  velocity  in  gas  particle  flow,  fps 

Up  = particle  velocity,  fps 

Up'  — nondimensional  particle  velocity 

11'  — mass  flow  rate  of  gas  particle  mixture  or  mass  flow  rate 

of  gas-only  flow  (ib/scc) 
x,y  — Cartesian  coordinates 

Si  l bscrijtts 

0 - pas 

P 3-.  particle 

1 — condition  at  starling  I'oinl  of  the  stream  tubes 


Introduction 

HP  II E study  of  gas  particle  flow  over  compressor  or  turbine 
blades  is  an  area  of  recent  interest.  This  is,  in  part, 
due  to  the  development  of  high  energy  propellants  for  com- 
bustion system  in  which  the  product  of  combustion  may  par- 
tially consist  of  finely  divide  i particles.  For  example,  tur- 
bines which  arc  invariably  present  in  nuclear  and  liquid  pro- 
pellant rockets,  are  subjected  to  gas  particle  flow.  Air- 
breathing  engines  operating  in  desert  areas,  where  the  inlet 
airflow  may  contain  sand  particles,  provide  another  example. 

The  flow  of  gas  particle-suspension1-*  involves  gas  particle 
interaction  through  viscous  drag  and  heat  transfer.  This 
interaction  causes  changes  in  the  properties  of  gas  and  par- 
ticles. These  particles  are  swept  through  the  cascade  nozzle 
by'  gas  flow,  lagging  behind  the  gas  in  velocity  and  tempera- 
ture. Both  effects  lead  to  the  deterioration  in  the  perform- 
ance of  compressor  or  turbine.  Furthermore,  the  particles 
may  cause  blade  erosion  and  consequently  failure.  Even 
when  the  erosion  is  not  severe  enough  to  cause  rupture, 
changes  of  blade  geometry  affect  cascade  performance. 

Blade  erosion  depends  upon  several  factors,  some  of  which 
are:  the  properties  of  the  blade  material;  the  total  mass  of 
impinging  particles;  the  particle  speed  and  the  angle  of 
attack;  and  the  temperature  of  both  the  gas  and  the  particles 
at  the  blade  surface.  Before  approaching  these  problems 
it  is  necessary  to  determine  the  way  in  which  the  gas  particle 
properties  vary  with  the  distance  along  the  surface  of  the 
blade. 

Method  of  Solution 

A numerical  solution  to  the  problem  of  gas  particle  flow 
over  blades  in  cascade  is  obtained,  utilizing  the  known  solu- 
tion for  isentropic  flow  without  particles  over  the  same  blade. 
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Pressure  Distribution  on  Blades  in 
Cascade  Nozzle  for  Particulate  Flow 

W.  Tabakoh  * and  M.  F.  llussuixt 
University  of  Cincinnati,  Cincinnati,  Oliio 

Nomenclature 

A «-  cross-scctionfil  area  of  stream  tube  (ft?; 

a *=  particle  concentration 

C,  = specific  heat  (Btu/lb’lt) 

dv  = particle  mean  diameter  (ft) 

g * gravity  constant 

J « mechanic  d equivalent  of  heat 

- coefficient  of  conductivity  for  the  g:es  (lltu/hr  ft°R) 

1*  “ gas  particle  suspension  pressure 

7>  = pressure  of  gas-only  flow 

Pr  « Frandtl  number 

l!j  *=  g:w$  constant 

A’tj  = Reynolds  number 

p — the  gas-only  flow  density  (lb/ft*) 

Pi?  — the  gas  density  in  g.u-;  particle  flow  (lb/ft*) 

Pp  -■*  the  particle  density  (lb/ft*) 

pr  — solid  particle  material  d*  usity  (lb/ft5) 

T =*  the  gas-only  flow  temperature  (°U) 

Tg  = the  gas  temperature  in  gas  par’icle  flow  (°R) 

u = gas-only  flow  velocity  (fps) 

v0  «=  gas  velocity  in  gas  particle  flow  (fps) 

= particle  velocity  (fps) 

W *=  mass  flow  rate  of  gas  particle  mixture  or  mass  flow  rate  of 
gas-only  flow  in  a stream  tube  (lb/sec) 

Superscripts 

' *=-  nondimension nl  parameter 

Subscripts 

g **  g-w  # 

p *=  particle 

l *-■  condition  at  starting  point  of  the  stream  lube 
0 = total  conditions 

Introduction 

r|'vIIE  pressure  distribution  on  the  blade  in  cascade  increases 
with  (he  introduction  of  solid  particles  to  the  gas  flow.1-- 
The  forces  on  the  blade  are  larger  thus  requiring  a stronger 
blade  design.  In  this  work,  a theoretical  approach  to  esti- 
mate the  pressure  distribution  on  the  b-adc  surface  for  a par- 
ticulate ges  flow  for  incompressible  and  compressible  cases 
was  developed.  Gas  flow  without  particles  past  a cascade 
was  considered  first.  It  was  assumed  that  an  experimental 
pressure  distribution  on  the  blade  surface  is  known  fur  given 
gas  flow  and  inlet  conditions.  It  was  assumed  that  two 
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stream  tubes  exist  in  the  flowfield  around  the  blade;  one  at 
the  suction  side  and  one  at  the  pressure  side.  The  pas  flow 
without  particles  was  used  to  determine  the  nundimensional 
area  of  the  stream  tube  as  function  of  the  given  pressure  dis- 
tributi 'ii  and  inlet  gas  conditions  of  the  nonp^rliculate  gas 
flow.  The  governing  equations  for  the  particulate  flow  were 
formulated  for  incompressible  and  compressible  eases.  The 
resulting  differential  equations  were  solved  numerically  for  the 
pressure  distribution  and  the  particulate  flow  properties. 

The  different  parameters  in  the  governing  equations  were 
nondiincnsionaJizcd  with  respect  to  values  at  a starting  point 
1,  which  is  in  the  vicinity  of  the  blade  loading  edge.  'J  hose 
parameters  are  expressed  as  follows  in  Kq.  (1) 

V ==  r/Ti  p'  = v/p  1 u'  « u/n,  p'  - p/p,  P'  = P/Pi 

U g — Ug/\tgl  7 , = 1 ii / f j,  Pa  Pa/ Ppt  U r tLp/tlgx 

p'r  - p,/p„  T'p  - Tp/T„  (1) 

A'  = A/A\  — 1/m'  (Incompressible  case)  (2) 

A ' A/Ai  *»  1/m'p'  (Compressible  case)  (C) 

’j'!ie  definition  of  the  particle  concentration  a is  expressed  ns 

a = ir,/(ir,  -i-  ir.)  - uyiv  (-D 

Incompressible  Mow 

The  continuity  and  momentum  equations  for  gas-only  flow- 
are 

(l/u)dn/ds  -r  (l/d}d.t/ds  - 0 (5) 

pudu/ds  = — (ydp/ds)  (6) 

Substituting  Eq.  (o)  into  Eq.  (0)  and  putting  the  resulting 
equation  in  nondimensiona!  form  by  using  Eq.  (1),  one  gets 
an  expression  for  the  rate  of  change  of  stream  tube  nondimen- 
siona) cross-sectional  area  as  a function  of  the  rate  of  change 
of  the  nondimcnsional  gas  without  particle?  flow  pressure. 

dA'/ds  = (A  '/u'!)  (ffp,/ pui*)  (dp'/  ds)  (7) 

where  u'  is  given  by 

u'  = [(;>„  - v/V o _ Pi)]1'5  (S) 

'i'be  equation  for  the  drag  force  on  the  solid  particles,  the 
momentum  equation  of  gas  and  particles,  and  the  continuity 
equation  of  gas  are 

UpdUp/ds  - 18;i„/ppdp!(»j  — uf)  (9) 

rf/’/ds  = — (ir/ff.l)  1(1  — a)du,/ds  + adUp/ds]  (10) 

p,utA  — - tl  — n)lV  (11) 

Differentiation  of  Eq.  (11)  and  substitution  of  l’qs.  (1)  and 
(7)  into  the  resulting  differential  equation  results  in 

du'./ds  = — (u',/.l')(d.l'/ds)  = 

- (u  V u'1)  (g/a/pui5)  (dp  '/ds)  (12) 

Substituting  Eq.  (1)  into  Kq,-'.  (0)  and  (10),  they  take  the 
nondimensiona!  form 

u'pdu'p/ds  = (18 fi,/u,xJ5pd,i){u',  - u'p)  (13) 
dV'/ds  ~ — p,u,*/gPiA'[du',/ds  -(■  (a/1  — cc)du',/ds]  (14) 

ii  4 


W.  Tabakoff  Dr.,  Professor,  Department  of  Aerospace  Engineering  Univ.  of 

Cincinnati,  Cincinnati,  Ohio,  U.S.A.  Memh.  A SME 
M.  F . Hussein  Graduate  Research  Asst.,  Department  of  Aerospace  Engineering, 
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The  gas-particle  flow  properties,  namely  flow  pressure,  qas  and  particles  velocities  and 
temperatures  were  calculated  theoretically  in  a b'ade-to-blade  stream  surface  of  a cascade 
nozzle.  The  effects  of  the  particle  concentration,  mean  diameter  and  materia!  density  on 
the  flow  properties  at  any  point  in  the  cascade  nozzle  were  investigated.  Some 
experimental  results  were  obtained  from  cascades  wind  tunnel  equipped  with  a solid 
particle  injecting  system.  The  increase  in  pressure  due  to  the  introduction  of  solid 
particles  gives  rise  to  forces  on  the  blades  indicating  that  this  information  is  of 
importance  to  the  turbornachinery  designer. 
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Hussein , M.F.  and  Tabakoff,  W.  , "Calculation  of  Particle  Tra- 
jectories in  a Stationary  Two  Dimensional  Cascade, " Project 
Themis  Report  72-27,  June  1972 . ( AD  7 G 4 2 6 7) 

Abstract 

The  trajectories  of  solid  particles  entrained  by  fluid  flow 
in  turbine  and  compressor  staters  were  calculated.  The  impact 
and  rebound  phenomena  of  the  solid  particles  was  investigated 
experimentally  and  then  considered  in  the  solution  of  the  equations 
of  motion  of  the  solid  particles.  The  effect  of  the  particle  mean 
diameter,  material  density,  and  initial  particle  and  gas  velocities 
on  the  dynamic  behavior  of  the  solid  particles  through  the  cascade 
channel  was  also  investigated.  In  addition,  this  study  yielded 
information  concerning  blade  erosion  damage. 


Hu  use  in , M.F.  and  Tabak  off,  V. . , "Calculation  of  the  Three  Pimen- 
si o n a 1_  Par t icle  Trajecto ri e s in  a Te r bine  Stage,"  Project  'i'hemis 
Report  No.  72-33,  September  1972.  (AD  7533G41 

Abstract 

The  equations  of  motion,  in  three  dimensions,  of  solid  particles 
entrained  by  a gas  flow  through  the  stationary  and  rotating  cascades 
of  a turbine  are  derived.  The  gas  velocity  components  and  density 
at  all  the  mesh  points  of  a square  grid  constructed  in  the  flow 
channels  are  computed  assuming  a compressible  flow.  Formulas  to 
determine  the  proper  drag  on  the  particles  for  a wide  range  of 
Reynolds  numbers  are  given.  A gas  particle  flow  tunnel  is  used  to 
investigate  experimentally  the  phenomenon  of  particle  impact  with 
the  turbine  blades  or  casing  and  their  rebound  from  these  walls. 
Formulas  for  the  restitution  ratio  due  to  collision  and  the  rebound 
to  incidence  angle  ratio  are  derived.  This  information  is  used  in 
the  equations  of  motion  of  the  solid  particles. 

The  dynamic  behavior  of  the  solid  particles  in  the  turbine 
stage,  namely  their  absolute  and  relative  trajectories,  absolute 
nondimensional  velocity  history  in  the  channel,  and  their  velocity 
diagrams  as  compared  to  that  of  the  gas,  is  investigated . The 
effect  of  different  flow  parameters,  mainly,  the  particle  mean 
diameter,  material  density,  and  particle  and  gas  initial  velocities 
on  the  dynamic  characteristics  of  the  solid  particles  are  studied. 

Observations  concerning  the  erosion  damage  suffered  by  the 
turbine  stator  and  rotor  blades  as  well  as  the  turbine  casing  due 
to  the  solid  particle  impingements  are  presented. 
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Abstract 

The  boundary  layer  of  gas-particle 
flows  in  fields  of  changing  pressures  is 
analyzed.  The  analysis,  which  utilizes 
the  momentum  integral  method,  can  be  used 
to  study  the  boundary  layer  of  any  partic- 
ulate flow  with  an  incompressible  gas- 
phase,  provided  the  flow  properties  out- 
side the  boundary  layer  are  known. 

Numerical  integration  of  the  two  momentum 
equations  was  carried  out  for  the  case  of 
particulate  flow  in  cascades.  For  this 
example,  the  results  indicate  that' the 
introduction  of  particles  leads  to  an  in- 
crease in  the  gas  boundary  layer  thickness. 
In  addition,  it  was  found  that  the  gas 
boundary  layer  characteristics  are  more 
sensitive  to  particle  concentration  than 
any  other  pe%rticulato  flow  parameter. 
However,  the  particle  boundary  layer 
characteristics  are  more  sensitive  to  mean 
particle  diameter  than  to  particle  con- 
centration. Other  parameters  such  as  inlet 
particle  gas  velocity  ratio  and  particle 
solid  material  density  have  more  pronounced 
effects  on  the  particle  boundary  layer 
than  on  that  of  the  gas. 
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This  paper  describes  the  results  from  an  investigation  of  the  gas-particle  flow 
trajectories,  velocities  and  pressure  distribution  in  an  axial  flow  turbine  stage. 
A gas-particle  flow  cascade  tunnel  and  high-speed  photographic  techniques  were 
used  to  conduct  the  experimental  investigation.  The  pressure  distribution  on  the 
blade  surface  was  measured  and  compared  with  the  theoretical  analysis,  the 
results  exhibiting  good  agreement  between  the  developed  theory  and  experiment. 
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IV'idt'n  Tabukoff')  and  Awatef  flamed  **),  Cincinnati,  Ohio/UOA  L)K  53. 529.5  : 532.526 

The  Boundary  Layc-  of  Pariicufats  Gas  Flov/+) 

fjfv  :■ :J:t : Die  Arbeit  befaft  sich  tr.it  drier  theorclitchen  Unlcrsuchuttg  tier  Inminnren  Gr eiizcchidit  fiir  ilic  vbenc 
Ga?-Partikel-C.  triinung.  Die  hr.pu'tintcgralglcichuiigcii  h i inkompressiblcr  Gasp’ase  und  bci  Vcrhandcmein 

c. ne$  Druckgradicnten  warden  nbgclcitct.  is  wird  dcr  full  drier  Mischstromung  ven  Gr.s  und  Partikelu  l ring'  cuter 
cbenen  Pl.itte  betrnhiet,  wobci  beicie  Pharen  dieselbe  Geschwindigkcit  in  der  Aupen>tromuug  hnben.  Die  tins 
einer  numeric  J:en  t cred'.nung  gewennenen  Ergcrnitse  iverden  in  Kurven  dnrgcstcllt,  in  denen  die  Crcnczmidit- 

d. cke  d s Gates,  die  Schubfpunn ung  an  der  IVnrid,  dcr  Formparumeter  fiir  this  Gcschwindigkcitsprofil  dc.  G.  es 
und  die  Vcrdr’ingurtgsdidce  d-  t Pat*'kcl  wiedergegt'en  sind.  Allc  diese  Parameter  sind  jeweils  fiir  vcrschiedcne 
Partikelkonzeui:  atiitnen  in  der  Anfvn-trdmung,  x'crschicdenc  Werte  dcr  Pc.rtikcldurchrnesser  und  vcrschiedcne 
IVcrte  dcr  Pi  Jit  e (let  fasten  TcilJten  dargestellt. 

Summary:  A theontic.il  analysis  of  the  laminar  boundary  layer  for  the  two-dimensional  gac-  particle  flew  is 
presented.  The  momentum  integral  etjuatuni:-  of  (Up  boundary  layer  with  incompressible  g».  phase  and  in  the 
presence  of  pr>  urc  gradient  are  derived.  The  uisc  of  a gas  particle  mixture  flow  oxer  a fit.:  plate  with  both 
phiv-es  hazing  the  same  velocity  in  th j free  stream  is  presented.  A numerical  calculation  is  carried  cut  and  the 
results  arc  presented  in  curi  es  representing  the  gas  boundary  layer  thickness,  the  shear  stress  at  the  wall,  Pc  gas 
velocity  profile  shape  factor , and  the  pevticlc  du  placement  thickness.  All  these  parameters  arc  presented  for 
various  particle  concentrations  in  the  free  sit  cam,  different  values  of  particle  diameter,  and  different  values  of 
solid  particle  material  densities. 

Resume:  L'itudc  perte  sur  u\e  analyse  theorique  tic  la  couehe  limite  lamir.aire  pour  I'ecoulcment  gaz-particule s 
bidimcnsionnel.  On  deduit  les  etp.u.tions  de  lii.tegrale  d'impulrion  tie  la  couche  limite  cn  cas  d'une  phase  de  gaz 
incompressible  et  en  presence  dun  gradient  dr  pression.  On  presents  le  cas  d'un  ccoulrrncnt  mis/e  se  composant 
de  gaz  et  dc  particules  le  long  d'une  plaque  plane,  les  deux  phases  ayani  la  mane  vitesse  dans  I'ecoulcment  libtc. 
Les  result. its  obtenus  d'un  calcul  nurnenque  sent  dcniontres  e:i  courbcs  rcpriscutant  Vepaisscur  de  couche  limite 
du  gaz,  la  tension  de  cisaillement  a la  paroi,  le  factcur  de  forme  pour  le  profit  dc  vitesse  du  gaz  cl  Vepaisscur  des 
particules.  Tons  ces  parametres  sent  detc  rrdnes  pour  diirerses  concentrations  des  paitic'dcs  dans  I'ecoulcment 
libre,  pour  diverges  valours  des  diatnclres  des  particules  et  pour  diverges  valeurs  de  la  densite  des  partieitlcs 
solides. 

1.  Introduction 

Cas  solid  flow  systems  are  of  considerable  importance  in 
the  present  rapidly  exp  mding  technology.  Typical 
examples  arc  the  two-pliase  flow  in  thrust  nozzles  using 
propellants  with  metal  additives,  pneumatic  conveying 
systems,  nuclear  reactors  with  gas  solid  feeding,  two- 
phase  flow  due  to  nbi  .lion  ceding,  and  many  other  .. 

The  concepts  of  continuum  media  ics  can  be  extended 
to  particulate  flow.  The  effect  of  a solid  boundary  on  a 
gas  particle  mixture  flowing  parallel  to  it  is  confined  to 
the  thin  region  near  the  wall.  The  concept  of  a viscous 
flow  boundary  layer  can  be  extended  to  a gas  particle 
flow.  Although  most  boundary  layers  with  particle 
suspension  are  turbulent,  the  case  of  a laminar  boundary 
layer  of  a gas  solid  suspension  is  studied  with  the  pur- 
pose of  developing  a basic  understanding  of  the  inter- 
action of  gas  solid  flow  with  the  wall.  Due  to  the 
complexity  of  the  boundary  layers  of  particulate  flows, 

•)  Frofcssor,  Department  of  Aerospace  Engineerinp,  University  of  Cincinnati. 

•*)  Graduate  Research  Assistant,  Department  of  Aerospace  Engineering,  University  of  Cincinnati, 
t)  This  research  work  was  supported  by  the  U.  5.  Army  Research  Office,  Durham,  under  Project  Themis  Contract  Number  DAHC  C4  -6,3-C-0Cl6. 


some  of  the  studies  done  on  such  problems  [1]  obtained 
the  first  few  terms  of  asymptotic  solutions  as  U/F  x->0 
which  corresponds  to  very  small  particle  sizes  or  very 
large  distances  from  the  leading  edge,  while  others  [2] 
obtained  also  the  first  few  terms  of  the  asymptotic 
solution  as  WF  x-*oc  which  corresponds  to  \v:y  large 
particle  sizes  or  very  short  distances  from  the  leading 
edge.  In  other  cases  [3]  the  gas  phase  boon  'ary  layer 
characteristics  were  assumed  to  be  unaffected  by  the 
presence  of  the  particle  to  obtain  an  estimate  of  the 
particle  phase  boundary  layer  thickness  very  close  to 
the  leading  edge.  Similar  approaches  were  used  for  un- 
steady problems  of  infinite  rotating  d'sk  [4]  and  both 
impulse  and  oscillating  motions  of  infinite  flat  plate 
[5,  6],  In  the  present  analysis  the  momentum  integral 
method  known  for  the  gas  only  flow  is  used  for  the  gas 
particle  flows.  The  two  resulting  ordinary  differential 
equations  are  coupled.  Here  the  momentum  integral 
equations  for  the  flat  plate  were  simultaneously  inte- 
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Ai  .STRACT 


A theoretical  method  was  developed  for  predicting  the 
pressure  distribution  over  a blade  in  cascade  for  a compres- 
sible flow  with  solid  particles.  Experimental  results  were 
obtained  ^rom  a ca5;cade  wind  tunnel  equipped  with  a solid 
particle  injection  system.  Good  agreement  v/as  noted  between 
the  theoretical  and  experimental  pressure  distribution.  The 
change  in  pressure  due  to  the  particles  gives  reduction  in 
the  force  on  the  biades. 

The  results  of  the  experimental  trajectories  investiga- 
tion can  provide  a good  approximate  prediction  of  the 
trajectories  and  velocities  of  solid  particles  suspended  in 
a fluid  and  passing  through  an  axial  flow  turbine  stage.  The 
technique  presented  allows  the  determination  of  the  blade 
areas  which  are  subjected  to  erosion  damage. 

The  presence  of  solid  particles  in  the  rocket  turbine 
gas  flow  changes  the  turbine  performance.  The  overall 
turbine  efficiency  decreases  as  a result  of  the  introduction 
of  solid  particles.  The  performance  experiment  was  performed 
on  an  oxidizer  pump  drive  turbine  for  an  M-l  rocket  engine. 
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s Tonal  Particle  Trajectories  in  a One  and  One-Half  Stage  of  a 
Compressor,"  Project  Themis  Report  No.  73-34,  February  1973. 

TAD  797164) 

Abstract 

The  equations  of  motion,  in  three  dimensions,  for  solid 
particles  moving  in  a compressible  gas  flow  in  a rotating  cascade 
are  solved  for  the  case  of  the  particles  moving  through  a com- 
pressor guide  vane  followed  by  a compressor  stage.  The  solution 
considers  the  particle  impact  and  rebound  with  both  the  blade  walls 
and  the  compressor  casing.  The  three  dimensional  absolute  paths 
of  the  particles,  their  trajectories  relative  to  the  compressor 
rotor,  their  velocity  histories,  and  the  combined  particle  velocity 
diagrams  in  the  compressor  cascades  are  given.  The  effects  of 
flow  parameters  on  the  dynamic  behavior  of  the  solid  particles 
throughout  the  compressor's  one  and  one-half  stage  are  investi- 
gated. Parameters  considered  included  particle  mean  diameter 
and  material  density,  initial  particle  and  gas  velocities  at  the 
guide  vane  inlet,  and  compressor  rotor  rotational  speed.  Obser- 
vations concerning  the  erosion  damage  to  the  blades  of  the 
compressor  cascades  are  also  discussed. 
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Abstract 


The  equations  that,  govern  the  three 
dimensional  motion  of  solid  particles  sus- 
pended by  a compressible  gas  flow  through 
a rotating  cascade  of  a turbomachine  are 
formulated.  These  equations  are  solved 
for  the  case  of  flow  through  a turbine 
stage.  The  solution  takes  into  account 
the  loss  in  particle  momentum  due  to  their 
collision  with  the  turbine  blades  or  cas- 
ing. The  dynamic  characteristics  of  the 
solid  particles;  namely,  their  absolute 
trajectories,  paths  relative  to  the  tur- 
bine rotor,  velocity  distributions,  and 
the  combined  stage  Velocity  diagrams,  are 
calculated.  The  effects  of  changing  the 
particles  mean  diameter,  material  density, 
and  initial  particle  and  gas  velocities  at 
the  stator  inlet  on  the  dynamic  character- 
istics of  the  solid  particles  are  investi- 
gated. The  results  obtained  from  this 
study  indicate  the  locations  on  the  tur- 
bine blades  subjected  to  severe  erosion 
d arr.a  ge . 
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ANALYSIS  OF  NONEQUILIBRIUM 
by 

A.  HAMED  and  W.  TAEAKOFF 
University  of  Cincinnati 
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ARTICULATE  FLOW 


Abstract 

The  equations  of  conservation  of  an- 
gular momentum  of  the  two  phases  of  a par- 
ticulate flow  were  derived,  and  the  stress 
tensor  was  determined.  It  was  found  to 
consist  of  anti-symmetric  elements  in  ad- 
dition to  the  pressure  and  symmetric  shear 
stresses  in  nonparticulate  flows.  An  im- 
plicit second  order  finite  difference 
scheme  was  used  to  determine  numerically 
the  particulate  flew  properties  throughout 
the  flow  field  from  frozen  to  equilibrium 
regimes.  The  numerical  study  investigated 
the  impulsive  motion  of  an  infinite  flat 
plate  in  an  incompressible  viscous  gas 
with  suspended  solid  particles.  The  ef- 
fect of  some  parameters  on  the  resulting 
particulate  flow  was  obtained. 
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Dynamic  Behavior  of  Solid  Particles  Suspended 
hy  Polluted  Flow  in  a Turbine  Stage 

M.  Fat  by  Hussein*  and  W.  Tabakofft 

University  of  Cincinnati,  Cincinnati,  Ohio 

The  equations  that  govern  tin*  three  dimensional  motion  of  solid  particles  suspended  by  a compress- 
ible gas  flow  through  a rotating  cascade  of  a turbomachine  are  formulated.  These  equations  are  solved 
for  the  case  of  flow  through  a turbine  stage.  The  solution  takes  into  account  the  loss  in  particle  mo- 
mentum due  to  their  collision  with  the  turbine  blades  or  casing.  The  dynamic  characteristics  of  the 
solid  particles;  namely,  their  absolute  trajectories,  paths  relative  to  the  turbine  rotor,  velocity  distri- 
butions, and  the  combined  stage  velocity  diagrams,  are  calculated.  The  effects  of  changing  the  par- 
ticles mean  diameter,  material  density,  and  initial  particle  and  gas  velocities  at  the  stator  inlet  on  the 
dynamic  characteristics  of  the  solid  particles  a^e  investigated.  The  results  obtained  from  this  study 
indicate  the  locations  on  the  turbine  blades  subjected  to  severe  erosion  damage. 


Nomenclature 

e = particle  absolute  acceleration 

li  = frame  fixed  in  blades 

$ = angle  between  particle  relative  velocity  and  tangent  to 

surface 

C = absolute  velocity 

C ~ absolute  velocity  component  in  the  x.6  plane 
Cpg  — specific  heat  of  gas  at  con  tant  pressure 
4'  - modified  stream  function 

Z>  = drag  force  on  spherical  part  icles 
dp  = particle  mean  diameter 

G = coefficient  inversely  proportional  with,  the  particle  char- 
acteristic time 

g(/.Y)  = Reynolds  number  dependent  function 
h - norma!  stream  channel  thicknc-s  (blade  height) 

X = angle  between  meridional  streamline  and  engine  axis 
nip  = mass  of  one  particle 
us  = gas  viscosity 

p ~ pressure  at  a point 

f*  = force  on  a spherical  parucic  due  to  pressure  gradient 

R = blade  mean  radius 

Re  - Reynolds  number 

r - radius  from  axis  of  rotation 

Pi  = gas  density 

i)it  = particle  material  density 

s = blade  angular  spacing 

T = temperature 

C = blade  speed  at  the  mean  radius 

u = relative  velocity  component  in  the  x-direciion 

V = relative  velocity 

V - relative  velocity  component  in  the  x.O  plane 

i = relative  velocity  component  in  the  tangential  direction 

W = weight  flow  rate  of  mixture  per  channel 
ic  = relative  velocity  component  in  the  radial  direction 
u?  = blade  angular  velocity 

x.P.r  = coordinate  axes  in  the  meridional,  tangential  and  radial 
directions,  fixed  in  B or  components  of  the  relative 
particle  position  vector 

X.  Y.Z  = axes  fixed  in  engine  or  components  of  the  absolute  parti- 
cle position  vector 
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x.b.z  — particle  velocity  components  measured  in  B 
x,  0,  2 = particle  acceleration  components  measured  in  B 
X.  Y.Z  = particle  absolute  velocity  measured  in  K 
X,  Y.Z  = particle  absolute  acceleration  measured  in  F. 
y - distance  measured  along  the  tangential  direction 

Subscripts 

a = absolute  trajectory 

g = gas 

I = initial  conditions  at  particles  entrance 

in  = conditions  at  the  boundary  A H 

n = normal  to  blade  surface 

p ~ particle 

t - tangent  to  blade  surface 

1 = before  collision 

2 — after  collision 


Introd  net  ion 

ROCKETS,  aircraft  engines,  and  industrial  gas  turbines 
operating  in  desert  areas  and  in  places  where  the  atmo- 
sphere is  polluted  by  small  solid  particles  can  be  examples 
of  machines  operating  under  gas  part icle  two  phase  flow 
conditions.  The  solid  particles  mixed  with  the  inlet  air  or 
combustion  gases,  due  to  the  difference  in  their  inertia,  will 
be  driven  away  from  the  streamlines  of  the  gas  and  impact 
with  the  surrounding  walls  of  the  engine.  The  presence  of 
solid  particles  in  the  Row  constitutes  a cause  for  severe  ero- 
sion damage  especially  to  the  rotating  parts  of  the  engine, 
where  the  particle  velocities  and  frequency  of  impacts,  as 
well  as  the  Row  temperature,  are  higher.  Hence,  turbines, 
axial  or  centrifugal,  of  rockets  or  gas  turbines  are  parts  of 
the  engine  critically  subjected  to  particles  erosion.  The 
centrifugal  forces  acting  on  the  particles  tend  to  force  them 
to  move  radially,  and  hence,  impact  with  the  turbine  casing, 
causing  it  to  also  suffer  from  erosion  damage.  In  order  to 
understand  and  predict  the  erosion  phenomenon  of  rotating 
turbomachines,  it  is  important  to  study  the  dynamic  char- 
acteristics or  behavior  of  the  solid  particles  entrained  by 
the  gas  flow  through  the  stages  of  turbines.  By  the  dynamic 
behavior,  it  is  meant,  the  absolute  and  relative  position  of 
the  particles  everywhere  in  the  channels,  their  velocity 
history  and  velocity  diagrams  as  well  as  a description  of  the 
collision  and  rebound  mechanisms  of  the  solid  particles 
from  the  walls.  The  effect  of  different  Row  parameters  such 
as  particle  mean  diameter,  material  density  and  initial  par- 
ticle and  gas  velocities  on  the  particle  dynamic  behavior  are 
investigated.  Results  of  this  study  are  then  used  to  make  ob- 
servations concerning  the  areas  of  the  blades  or  casing  that 
are  subjected  to  erosion  damage. 

The  three  dimensional  equations  of  motion  of  solid  parti- 
cles moving  in  a compressible  gas  stream  in  a rotating  cas- 
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Hamed,  A.  and  Tabakoff , W. , "The  Solution  of  the  Equations  of 
Motion  for  Unsteady  Viscous  Particulate  Flows,"  Report  No. 

74-42,  July  1974. 

Abstract 

A numerical  method  is  developed  to  solve  the  equations  that 
govern  the  motion  of  the  unsteady  two  dimensional  flow  of  a sus- 
pension consisting  of  incompressible  viscous  gas  and  suspended 
solid  particles.  The  numerical  method  is  implicit,  second  order 
and  iterative.  The  frozen  and  equilibrium  wall  friction  co- 
efficients are  compared  with  the  corresponding  analytical  deter- 
mined values  for  a simple  example,  and  are  found  to  be  in  agreement. 
It  is  concluded  that  the  present  numerical  method  for  investigating 
the  nonequilibrium  flow  regime,  gives  accurate  results  not  only 
for  nonequilibrium  flow  conditions,  but  also  predicts  frozen  and 
equilibrium  flow  conditions  very  accurately. 


Ball,  R.  and  Tabakoff,  W. , "An  Experimental  Investigation  of  the 
Particle  Dynamics  of  Quartz  Sand  Impacting  6A1-4V  Titanium  and 
410  Stainless  Steel  in  an  Erosive  Environment,"  Report  No.  74-43, 
October  1974. 

Abstract 

The  impact  and  rebound  characteristics  of  high  speed  erosive 
quartz  particles  have  been  experimentally  determined.  The  impact 
parameters  were  found  to  be  statistical  in  nature  and  the  sta- 
tistical distributions  were  obtained.  The  alloys  used  in  this 
investigation  were  410  stainless  steel  and  6A1-4V  titanium. 
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Some  Effects  Caused  by  Solid  Particles  in  Flows 

A.  HaMED*  AND  W.  TABAKOFFt 
University  of  Cincinnati,  Cincinnati,  Ohio 


Theme 

TPHE  equations  of  conservation  of  angular  momentum  of  the 
X two  phases  of  a particulate  flow  are  derived,  and  the  stress 
tensor  determined.  It  was  found  to  consist  of  antisymmetric 
elements  in  addition  to  the  pressure  and  symmetric  shear  stresses 
in  nonpar  liculate  flows.  The  particulate  flow  properties  are 
determined  numerically  throughout  the  flowfield  from  frozen  to 
equilibrium  regimes. 

Contents 

A study  of  the  behavior  of  particulate  flows  near  the  solid 
boundaries  is  necessary  to  determine  the  effect  of  the  presence  of 
solid  particles  on  the  skin  friction  and  heat  transfer.  In  the 
present  analysis,  the  particles  are  continuously  distributed 
throughout  the  flowfield  and  each  of  the  two  phases  has  its  own 
mean  properties.  The  particles  translational  velocity  is  generally 
different  in  magnitude  and  direction  from  that  of  the  gas.  They 
also  have  a rotational  velocity  which  is  not  equal  to  the  gas 
rotational  speed.  The  particles  motion,  however,  is  basically  due 
to  the  gas  motion. 

If  the  rate  of  change  of  the  angular  momentum  of  the 
suspension  in  the  volume  V is  set  equal  to  the  total  moment 
of  the  surface  forces  about  an  axis  in  the  k-diiection  passing 
through  a point  0 within  this  volume  we  obtain 

\0--£ipt.kllri(AvJdl)AV  +-J  /.prckllri(dvrJdt)dV  + 

J XP,r fidrujdt)  dV  = J E„>[d(r,<rIJ)/flr1]dl'/  (1) 

where  p and  pp  are  the  fluid  and  solid  particle  material  densities, 
■/  is  the  volume  occupied  by  the  solid  particles  per  unit  volume  of 
mixture,  ry  is  the  radius  of  gyration  of  a solid  particle,  i;  and  vp 
arethegasandparticlevelocities,cn  is  the  solid  particle  rotational 
velocity,  and  r is  the  position  vector,  whereas  the  subscripts  i,j,  k 
refer  to  components  in  these  directions. 

If  the  volume,  V,  is  reduced  to  zero  such  that  the  configuration, 
which  is  made  lip  of  the  boundary  of  the  volume  and  the  fixed 
point  0,  retains  the  same  shape,  the  last  term  on  the  left-hand 
side  and  the  first  term  on  the  right-hand  side  of  Eq.  ( 1 ),  approach 
zero  in  the  same  order  as  V,  while  the  rest  of  the  terms  approach 
zero  as  K4,J.  Thus 

XP,  r,  2( dwjdt ) = euj  c.j  (2) 

The  stress  tensor  is  written  below  ns  the  sum  of  three  terms. 
The  first  term  is  spherically  symmetric  and  represents  uniform 
compression,  the  second  is  a symmetrical  tensor,  and  the  third  is 
antisymmetric 

(Ty-  -p^+ty  + Ty*  (3) 
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where  p is  the  pressure  and 


If  the  expression  for  the  stress  tensor  given  by  Eq.  (3)  is  sub- 
stituted into  Eq.  (2),  the  contribution  of  the  spherically  symmetric 
and  symmetrical  tensors  to  the  right-hand  side  of  Eq.  (2)  is  zero 
and,  hence,  the  following  relation  is  obtained; 

r,/  = (xpflyAdioJdt) ; i * j ¥>  k (5) 

Thus,  the  stress  tensor  in  suspension  flows  consists  of  an  anti- 
symmetric shear  stress,  {t  f),  which  is  equal  to  half  the  rate  of 
change  of  the  local  angular  momentum  of  the  solid  particles  per 
unit  volume,  in  addition  to  the  familiar  pressure  and  symmetric 
shear  stress  in  nonparticulate  fluid  flows. 

Some  numerically  determined,  incompressible  particulate  flow 
properties  due  to  the  impulsive  motion  of  an  infinite  flat  plate1 
are  presented.  The  effect  of  the  initial  particle  concentration,  yc 
and  the  particle  to  gas  density  ratio,  pjp,  on  th-  flow  properties 
was  investigated.  The  friction  coefficient  is  plott  d against  time 
for  different  particle  to  gas  density  ratios  in  Fig.  1.  The  non- 
dimensional  skin-friction  and  antisymmetric  stress  tensor  are 
defined  as 


C/  = C/Re  t*)1/J 

r 2x„‘(RetT1 
pV2X* 

The  Reynolds  number.  Re.  is  based  on  the  plate  velocity,  U,  and 
the  relaxation  time  of  particle  translation,  r, 

Re  =-  U2rjv 

The  time  is  normalized  with  respect  to  r, 
t*  = 1/ r, 

where 

t,  = (d2/lSv)(pfp) 

In  the  last  equation,  d is  the  solid  particle  diameter  and  v is  the 
kinematic  viscosity. 


Fig.  t Effect  of  particle  to  gas  density  ratio  on  Ibe  friction  coefficient 

at  the  wall. 
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A NUMERICAL  METHOD  FOR  THE  SOLUTION  OF 


PARTICULATE  FLOW  EQUATIONS 
by 

A.  HAMED  and  W.  TABAKOFF 
University  of  Cincinnati 
Cincinnati,  Ohio 


Abstract 

A numerical  method  is  developed  to 
solve  the  equations  governing  the  motion 
of  a gas-solid  particle  suspension.  The 
values  of  the  skin  friction  coefficient 
at  the  plate  surface  at  zero  time  and  at 
time  equal  to  fifteen  times  the  charac- 
teristic time  of  particle  translational 
relaxation  are  compared  to  the  skin  fric- 
tion coefficients  of  the  frozen  and 
equilibrium  flow  regimes  determined  analyt- 
ically. The  two  values  arc  found  to  be  in 
agreement.  It  is  concluded  that  the 
present  numerical  method  for  investigating 
the  nonequilibrium  flow  regime,  gives  ac- 
curate results  not  only  for  nonequilibrium 
flow  conditions,  but  also  predicts  frozen 
and  equilibrium  flow  conditions  very  ac- 
curately . 
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COMPUTATION  AND  PLOTTING  OF  SOLID  PARTICLE 
FLOW  IN  ROTATING  CASCADES* 

M.  Fathy  Hussein  and  W.  Tadakoff 
University  of  Cincinnati,  Cincinnati,  Ohio  45221,  U.S.A. 


(Received  3 January,  revised  29  November  1973) 


Abstract — A numerical  technique  to  solve  the  three  dimensional  equations  of  motion  of  solid 
particle-;  suspended  by  a compressible  fluid  flow  in  a rotating  cascade  of  a turbomachine  is 
presented.  The  solution  of  these  equations  utilizes  the  numerically  computed  nonparticulatc 
fluid  flow  properties  through  the  same  cascade  channel.  It  also  takes  into  consideration  the 
experimentally  investigated  phenomenon  of  solid  particle  impingement  with  the  blade 
surface  and  turbomachinc  casing  and  their  subsequent  rebound.  Computer  plottings  of  the 
absolute  paths  of  the  solid  particles,  their  trajectories  relative  to  the  rotor  blades,  and  their 
velocity  distributions  in  a turbine  stator,  turbine  stage,  and  a compressor  one  and  one-half 
stage  arc  given. 


NOMENCLATURE 


P 

B 

C 

CD 

C" 

C,„ 

'K 

dr 

Ax 

A9 

A t 

E 

1 

G 

h 

A 

/».  _ 

Hi,  Hi,  "j 

R 

Rc 


P, 

Pt 

s 

T 

u 

V 

v 

W 
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angle  between  particle  relative  velocity  and  tangent  to  surface 
frame  fixed  in  blades 
absolute  velocity 

drag  coefficient  on  spherical  particles 
specific  heat  of  gas  at  constant  pressure 

particle  normal  velocity  component  relative  to  the  blade  at  impingement 
modified  stream  function 
particle  mean  diameter 

spacing  between  adjacent  points  in  the  axijl  direction 
spacing  between  adjacent  points  in  the  tangential  direction 
increament  of  time 
frame  fixed  in  the  engine 

angle  between  relative  velocity  and  meridional  plane 
coefficient 

normal  stream  channel  thickness  (blade  height) 
angle  between  meridional  streamline  and  engine  axis 
gas  viscosity 

unit  vectors  in  the  direction  of  the  coordinate  curves  x,  0,  z 

unit  vectors  in  the  direction  of  the  axes  X,  Y,  Z 

blade  mean  radius 

Reynolds  number 

radius  from  axis  of  rotation 

gas  density 

particle  material  density 
blade  angular  spacing 
temperature 

relative  velocity  component  in  the  .v-direction 
relative  velocity 

relative  velocity  component  in  the  tangential  direction 

weight  flow  rate  of  mixture  per  channel 

relative  velocity  component  in  the  radial  direction 
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Performance  and  Flow  Properties  Change  Through  a Rocket  Turbine  by 

Presence  of  Solid  Particles1 


Widen  Tabakoff',  Wishaa  M.  Hosny3  and  Awatcf  A.  I lamed4 
University  of  Cincinnati,  Ohio  45221,  U.S.A. 

(Received  25  September  1973) 

Abstract  Performance  and  Flow  Properties  Change  Through  a Rocket  Turbine  by  Presence  of  Solid  Particles. 
A theoretical  method  was  developed  for  predicting  the  pressure  distribution  over  a blade  in  cascade  for  a 
compressible  How  with  solid  particles.  Experimental  results  were  obtained  from  a cascade  wind  tunnel 
equipped  with  a solid  particle  injection  system.  Good  agreement  was  noted  between  the  theoretical  and 
experimental  pressure  distribution.  "I  he  change  in  pressure  due  to  the  particles  gives  reduction  in  the  force 
on  the  blades. 

The  results  of  the  experimental  trajectories  investigation  can  provide  a good  approximate  prediction  of  the 
trajectories  and  velocities  or  solid  particles  suspended  in  a fluid  and  passing  through  an  axial  flow  turbine 
stage.  The  technique  presented  allows  the  determination  of  the  blade  areas  which  are  subjected  to  erosion 
damage. 

The  presence  of  solid  particles  in  the  rocket  turbine  gas  flow  changes  the  turbine  performance.  The  overall 
turbine  efficiency  decreases  as  a result  of  the  introduction  of  solid  particles.  The  performance  experiment 
was  performed  on  an  oxidizer  pun.p  drive  turbine  for  an  M-l  rocket  engine. 

Translated  abstracts  appear  at  the  end  of  this  article 


Nomenclature 

A cross-sectional  area  of  stream  tube  (ft2); 

A'  non-dimensional  cross-sectional  area  of 
stream  tube; 

a the  ratio  of  the  mass  flow  rate  of  particles 

to  the  total  mass  flow  rate  of  the  gas  and 
particle  mixture; 

b.h  p.  brake  horsepower  (h.p.); 

Cpp  specific  heat  of  solid  particles  material 

(Btu/lb  CR); 

Cpg  gas  specific  heat  at  constant  pressure 

(Btu/lb  °R) ; 

<5s  distance  along  airfoil  contour  between 

two  successive  points  (ft); 

<5  ratio  of  inlet  total  pressure  to  NACA 

standard  S.L.  pressure  (P.,/14.696); 

dp  particle  mean  diameter  (ft); 


1 This  work  was  sponsored  by  the  U.S.  Army  Research 
Office - Durham,  under  Contract  Number  DAHC04-69C- 
0016. 
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3 Graduate  Research  Assistant,  Department  of  Aero- 
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4 Research  Assistant  Professor,  Department  of  Aero- 
space Engineering,  University  of  Cincinnati. 


g gravity  constant; 

O'  ratio  of  inlet  total  temperature  to  NACA 

standard  S.L.  temperature  (7Jt/519); 
h • enthalpy  (Btu/lb); 

J mechanical  equivalent  of  heat ; 

kg  coefficient  of  conductivity  for  the  gas 
(Btu/hr  ft  °R); 

N turbine  rotative  speed; 

Nn  turbine  design  speed ; 

P gas  particle  suspension  pressure; 

P'  non-dimensional  gas  particle  suspension 

pressure; 

p pressure  of  gas-only  flow; 

p'  non-dimensional  gas-only  flow  pressure; 

Pr  Prandtl  number; 

Rg  gas  constant; 

Re  Reynolds  number; 

P the  gas-only  flow  density  (lb/ft3); 

p'  the  gas-only  flow  non-dimensional  density : 

pg  the  gas  density  in  gas  particle  flow  (lb/ft3); 

p'g  the  non-dimensional  gas  density  in  gas 

particle  flow; 

pp  the  particle  density  (lb/ft3): 

p'p  the  non-dimensional  particle  density; 

Pp  solid  particle  material  density  (lb/ft3); 

T the  gas-only  flow  temperature  (°R)  or 
turbine  torque : 
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Solid  Particle  Demixing  in  a Suspension  Flow 
of  Viscous  Gas 

The  motion  of  a medium  consisting  of  an  incompressible  viscous j fluid  and  suspended 
solid  particles  was  analyzed  from  t/w  standpoint  of  continuum  mechanics.  It  was 
assumed  that  the  particles'  translational  and  rotational  velocities  were  different  fro  in 
those  of  the  fluid.  A numerical  scheme  was  used  to  determine  the  non-equilibrium 
particulate  flow  properties  as  well  as  their  equilibrium  and  frozen  values.  The  results 
arc  presented  for  the  case  of  particulate  flow  due  to  the  impulsive  motion  of  an  infinite 
flat  plate  in  a suspension.  .1  de  mixed  region , with  no  particles  present , was  found  to 
develop  near  the  plate  due  to  particle  migration  away  from  the  wall.  Similar  de  mixed 
particle  regions  were  noticed  in  the  experimental  data  of  particulate  flows  in  pipes 
and  channels. 


Introduction 


Systems  of  fliti'.l  flows  with  suspended  solid  particles  have 
practical  applications  in  numerous  industrial  areas.  In  some  of 
the  applications  the  contact  of  the  particles  with  the  solid  bound- 
ary should  be  avoided  in  order  to  reduce  erosion  damage,  as  in 
the  case  of  particulate  flows  in  nozzles  (l]1  and  turboriiachines 
[2],  or  the  transport  of  solid  material  in  suspended  flows  (3J.  In 
other  applications  like  filtration  and  separation  of  solid  particles, 
the  deposition  of  the  particles  on  the  wall  is  desired  [4).  Since  the 
effects  of  fluid  viscosity  also  prevail  near  the  boundaries,  the 
study  of  particulate  flows  with  a viscous  fluid  phase  is  necessary 
in  order  to  analyze  the  solid  particle  motion  at  the  walls. 

For  particulate  flow  systems  with  a nonviscous  fluid  phase, 
the  problem  is  basically  the  study  of  tiie  relaxation  of  transla- 
tional motion  and  temperature  of  the  particles  from  their  initial 
values,  to  the  fluid  velocities  and  temperature.  If  the  fluid  phase 
is  viscous,  the  relaxation  of  the  particle  rotation  as  a solid  body, 
which  is  initially  different  from  the  fluid  rotation,  should  be  in- 
cluded in  the  analysis.  A lift  force  of  interaction  exists  in  addi- 
tion to  the  drag. 

Although  many  of  the  fluid-solid  particle  suspension  flows  are 
turbulent,  theoretical  treatment  of  such  flows  is  formidable  since 
an  extensive  amount  of  experimental  data  is  needed  in  order  to 
obtain  semiempirical  hypothesis  necessary  to  treat  the  problem 


perimental  data  of  turbulent  particulate  flows  (lj,  [3],  [5],  |G] 
are  insufficient  to  arrive  at  a conclusive  hypothesis.  The  study 
of  laminar  suspension  flows  points  out  the  important  flow  param- 
eters arid  determines  their  effect  on  the  behavior  of  the  two 
phases. 


Analysis 


Governing  Equations.  The  equations  of  motion  are  written 
from  the  continuum  point  of  view,  for  a two  dimensional  fluid- 
solid  particle  medium,  treating  it  as  consisting  of  two  continuous 
homogeneous  phases,  each  having  its  mean  properties  [7],  (S), 
[9],  [10).  The  particles  arc  assumed  to  be  of  homogeneous  size 
and  physical  properties. 

The  equation  of  conservation  of  mass  of  the  fluid  is: 


dx 

dt 


31(1  - x)u] 


Old  - xM 


= o 


dx  , dy 

while  conservation  of  mass  of  the  particle  phase  is  given  by: 

dlxnj  , d{xr 

Ox 


d) 


+ 

at  ^ 


+ 


dy 


(2) 


The  equations  of  conservation  of  linear  momentum  of  the  fluid 
in  the  x and  y directions  arc: 


mathematically.  The  amount  and  accuracy  of  the  available  ex- 
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PARTICLES  IN  A CASCADE 

M.F.  HUSSEIN,  W.  TABAKOFF  and  A.  HAMED 
Department  of  Aerospace  Engineering 
University  of  Cincinnati 
Cincinnati,  Ohio  45221,  USA 

1.  Introduction 

This  investigation  was  undertaken  to  determine  a theo- 
retical approach  fo  calculate  the  dynamic  behavior  of  solid 
particles  entrained  by  a gas  flow  in  a two-dimensional 
stationary  cascade  of  a turbine. 

The  equations  of  motion  of  solid  particles  moving  in  a 
stream  of.  gas  were  formulated  assuming  -that  the  only  force 
exerted  on  the  particles  is  the  drag  force  which  causes  the 
particles  acceleration.  In  order  to  determine  the  solution 
of  the  equations  of  motion,  information  about  the  drag 
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Numerical  Method  for  Solution  of 
Particulate  Flow  Equations 


A.  Hamed*  and  W.  TabakoffI 
University  of  Cincinnati,  Cincinnati,  Ohio 


Theme 

THE  equations  describing  the  motion  of  a medium  consisting 
of  viscous  fluid  and  suspended  solid  particles  are  nonlinear, 
partial  differentia)  equations,  and  are  always  coupled  even  in  the 
case  of  incompressible  fluid  phase.  Analytical  solutions  could  be 
obtained  for  frozen  and  equilibrium  flow  regimes  of  a few 
problems  with  simplified  models  for  the  forces  of  bulk  interaction 
between  the  two  phases  If  the  nonequilibrium  particulate  flow 
is  to  be  studied,  numerical  computational  methods  have  to  be 
used.  A numerical  method  to  analyze  the  unsteady  two- 
dimensional  motion  of  incompressible  viscous  gas  and  suspended 
solid  particles  is  presented.  This  method  is  used  to  study  the 
particulate  flow  due  to  the  impulsive  motion  of  an  infinite  flat 
plate  in  an  otherw  ise  stationary  suspension. 


Contents 

The  motion  of  the  fluid  and  solid  particles  is  governed  by 
the  equations  of  continuity  and  momentum  of  the  gas,  the 
equations  of  linear  and  angular  momentum  of  particles,  and 
the  equation  of  conservation  of  mass  of  particles.  The  equation 
'of  angular  momentum  of  the  suspension  was  derived  by  the 
authors'  and  was  used  to  determine  the  stress  tensor  in  the 
particulate  flow  . The  equations  of  motion  are  solved  to  determine 
the  following  particulate  flow  properties;  the  fluid  velocity,  the 
particle  concentration,  and  the  translational  and  rotational 
speeds  of  the  solid  particles  throughout  the  llowficld.  The  solid 
particle  translational  and  rotational  velocities  are  generally 
different  from  those  of  the  gas  phase,  depending  on  the  flow 
regime.  In  the  frozen  or  near  frozen  regime,  the  difference 
between  the  translational  velocity  of  the  two  phases  which  is 
referred  to  as  the  slip  velocity  is  very  large.  The  particulate  flow 
is  said  to  be  in  the  near  equilibrium  and  equilibrium  regimes 
when  the  slip  velocity  is  small  and  approaching  zero.  The  non- 
equilibrium flow  regime  corresponds  to  the  intermediate  range 
of  slip  velocity. 

The  equations  of  motion  are  solved  in  nondimensiona!  form, 
the  normalized  suspension  flow  variables  are  defined  as 
t*  - tit,.  i i*  = u/U,  u‘  = up!U,  u>*  = [tI/(ffe),,:]u> 

= [fRe)Xlly! l/r,].  v*  = {Re)'2v/u;  vp*  = (Re)':1v„/U( 
where  Re  = U1  t,/v.  The  time  t is  normalized  with  respect  to 
the  characteristic  time  of  particle  translational  motion  r,  where 
t,  = (d2/l8v)lj>p,  (>).  u and  r are  the  gas  velocity  components  in 
the  v and  y directions,  up  and  vp  are  the  particle  velocity 
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Dynamics  Conference.  Palo  Alto,  California.  June  17-19,  1974;  sub- 
mitted June  IS.  1974;  synoptic  received  September  12.  1974;  revision 
received  November  11,  1974.  Full  report  on  which  Paper  74-561  was 
based  is  available  from  National  Technical  Information  Service, 
Springfield,  Va„  22151,  as  N'75-12227  at  the  standard  price  (available 
upon  request).  This  work  was  sponsored  under  Contract  DAHC04- 
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components  in  the  ,x  and  y directions,  ui  is  the  solid  particle 
rotational  velocity,  p the  gas  density,  pp  the  solid  particle 
material  density,  d the  solid  particle  diameter,  and  U is  a 
characteristic  velocity  of  the  flow. 

Before  solving  the  governing  equations,  a transformation  of 
the  independent  variables  is  used.  Self-similar  solutions  cannot 
be  obtained  in  the  case  of  particular  flows  due'  to  the 
characteristic  relaxation  times  associated  with  the  translational 
and  rotational  motions  of  the  particles.  A transformation  that 
leads  to  a self-similar  solution  in  nonparticulate  flows  is  used 
to  eliminate  the  discontinuity  that  would  otherwise  exist  in  the 
initial  gas  velocity  profile  in  the  physical  plane.  The  trans- 
formation also  makes  it  unnecessary  to  add  more  mesh  points 
in  the  direction  normal  to  the  flow  as  the  thickness  of  the 
momentum  boundary  layer  increases.  For  the  example  con- 
sidered here,  the  following  transformation  is  used: 

n = y*/ (2) 

The  forces  and  torque  of  bulk  interaction  between  the  two 
phases  depend  generally  on  the  range  of  some  parameters  of  the 
particulate  flow  under  consideration.  The  drag  force  and  torque 
due  to  the  difference  between  the  translational  and  rotational 
velocities  of  the  particles  and  the  fluid,  as  well  as  the  lift  force 
due  to  the  translational  motion  of  the  particles  in  the  shear  flow 
are  considered.  The  expressions  used  in  the  present  example 
correspond  to  dilute  suspensions  and  to  low  slip  Reynolds 
number,  ignoring  the  gravity  forces.  The  slip  Reynolds  number 
is  based  on  the  slip  velocity  between  the  particles  and  the  gas 
and  the  solid  particle  diameter.  The  numerical  method  presented 
here  can  be  used  however  for  any  other  range  of  flow  para- 
meters if  the  corresponding  expression  of  the  forces  and  torque 
of  bulk  interaction  are  used. 

The  equations  that  govern  the  unsteady  two-dimensional 
motion  of  the  suspension  in  the  transformed  plane  (>/.  r*|  are 

f[n -*>'*]  + c[y>*3  = 0 (3) 

crj  di/ 


2r*  ?-  + 2 (t 


*)'■  V cn  <>/  / 


T— + 

' v*  ii  \ du* 

[a* 

,2 (f*)1,T  ~ 2t*J 

1 c-V  1.5  d/du* 

— — [ 1 + 1.5/i]  -TT  + — ~t b 

4r*  L ay  4 1*  o;  di] 


3 c(ym«) 

2(t*)''2  ft; 


~%G(u*- 

P 


123  / \ Jveff'V  2 

2n  Z\  (Ret*)'  2 p dij ) 


+ (JA JL)  e<-  = G(u'-u  *)- 

Cl*  \2lr*)'  2 2 (*/  dii  p 

3»M'Y-y  i 

It  W \ -Hi*)' 2 dll ) 

iHff-W.f-— - 

It  w V 4U*)l2d,,. 


|w*-u  „*) 


Hussein,  M.F.  and  Tabakoff,  W. , "Computer  Program  for  Calculations 
of  Particle  Trajectories  Through  a Rotating  Cascade,"  Report 
No.  76-47,  May  1976. 

Abstract 

This  report  gives  a listing  of  the  computer  program  used  to 
calculate  the  three  dimensional  absolute  paths  of  particles 
through  a rotating  cascade,  their  trajectories  relative  to  the 
rotor  and  their  velocity  histories.  The  program  considers  the 
impact  and  rebound  of  particles  with  the  blade  walls  or  casing. 

The  program  is  written  for  axial  flew  machines  and  may  be 
modified  to  handle  other  types.  The  program  procedures,  descrip- 
tion of  the  program  input  and  output  as  well  as  program  sub- 
routines and  their  functions  are  discussed. 
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Effect  of  Solid  Particles  on  Turbine 
Performance 

A theoretical  method  was  developed  fur  predicting  the  pressure  distribution  over  a blade 
in  cascade  for  a compressible  flow  with  solid  particles.  Experimental  results  were  ob- 
tained from  a cascade  wind  tunnel  equipped  with  a solid  particle  injection  system  Good 
agreement  uas  noted  between  the  theoretical  and  experimental  pressure  distribution. 
The  change  in  pressure  due  to  the  particles  gives  reduction  in  the  force  on  the  blades. 
The  presence  of  solid  particles  in  air  breathing  engine  gas  flow  changes  the  turbine  per- 
formance. The  overall  turbine  efficiency  decreases  as  a result  of  the  introduction  of  solid 
particles  The  performance  experiment  uas  performed  on  a two-stage  velocity-com- 
pounded turbine. 


Introduction 

The  problem  of  a gas-flow  mixed  with  solid  particles  in  axial 
flow  turbines  has  great  importance  in  industrial,  naval,  and  aero- 
nautical applications.  For  example,  the  operation  of  engines  in 
desert  area*  or  solid  particle  polluted  atmospheres.  In  some  com- 
bustion processes,  solid  particles  are  formed  as  by-products  of  the 
combustion  This  is  especially  true  for  the  more  advanced,  power- 
ful fuels  The  presence  of  solid  particles  can  cause  erosion  damage 
of  the  turbine  blades  and  alter  the  pressure  distribution  on  them. 
The  erosion  damage  on  turbine  blades  from  a particulatcd  flow 
may  be  very  critical.  This  damage  is  especially  critical  in  turbine 
no?zles  where  the  temperature  of  the  gas-flow  and  the  solid  parti- 
cles are  high,  along  with  high  particle  velocities  which  greatly  in- 
crease the  erosion  damage.  Air  flow  with  suspended  sand  particles 
would  also  affect  the  compressor  and  the  turbine  performances 
and  consequently  the*  engine  performance. 

In  this  work  the  following  two  basic  problems  are  presented: 

1 A theoretical  analysis  which  is  used  to  estimate  the  pressure 
distribution  on  the  blades  for  a gas  particle  suspension  flowing 
through  a turbine  cascade. 

2 An  experimental  study  to  show  the  effect  of  particulatcd 
flow  on  a jet  engine  turbine  performance. 

Pressure  Distribution  on  the  Blades  for  Gas  Particle 
Flow 

In  references  (l,  2]1  a discussion  and  analysis  of  the  flow  proper- 
ties and  pressure  distribution  of  gas  particle  suspension  over  blade 


1 Numbers  in  brackets  designate  References  at  end  of  paper. 

Contributed  by  the  G a?  Turbine  Conference  Division  and  presented  at 
the  Gas  Turbine  Conference,  Houston.  Texas.  March  2-6.  1975,  of  THK 
AMERICAN  SOCIETY  OF  MECHANICAL  ENGINEERS.  Manuscript  re 
ceived  at  ASME  Headquarters  December  2.  1974.  Paper  No.  75-GT-41. 


surfaces  of  a cascade  nozzle  were  presented.  In  these  references  the 
gas  flow  was  considers-  ! inviscid  except  for  the  drag  it  exerts  on  the 
particles. 

The  theoretical  approach  to  the  solution  of  this  problem  is  to 
consider  gas  flow  without  particles  past  a cascade  to  determine  the 
pressure  distribution  over  the  blade  and  inlet  flow  conditions. 
Next,  it  was  assumed  that  two  stream  tubes  exist  in  the  flow  field 
around  the  blade;  one  at  the  suction  side  and  one  at  the  pressure 
side.  The  gas  flow  without  particles  was  used  to  determine  the 
nondimensional  area  of  the  stream  tube  ns  a function  of  the  given 
pressure  distribution  and  inlet  gas  conditions  of  the  nonparticulate 
gas  flow.  The  governing  equations  of  the  particulate  compressible 
gas  flow  were  formulated.  These  equations  were  solved  numerical- 
ly for  the  pressure  distribution.  The  different  parameters  in  the 
governing  equations  were  nondimensionalized  with  respect  to 
values  at  some  starting  point  which  is  in  the  vicinity  of  the  blade 
leading  edge. 

These  parameters  arc  respectively,  nondimensional  temperature 
(T'J,  pressure  (p’)t  velocity  (u’)t  and  density  ip')  of  the  gas  without 
particles  flow.  The  nondimensional  parameters  for  the  gas  particle 
flow  are:  particulate  flow  pressure  (P’)t  gas  velocity  (u\),  tempera- 
ture (T'K)  and  density  (f/K)  and  nondimensional  particle  velocity 
(u’p),  temperature  (V p)  and  density  (p'p)  for  gas  particle  suspen- 
sion. 


The  nondimensional  area  distribution  of  the  stream  tube  of  the 
gas-flow  on  the  airfoil  surface  can  be  calculated  from  the  continui- 
ty equation.  For  the  compressible  flow 
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Badnmv  dvnamiki  czastek  stalych  ; ich  crozyjne  dzialanie  w palisadzie 

Streszczenie 

Oblicza  si?  trajektorie  i pnjdkosci  cz.ystck  stalych  rtiesionych  przez  przeplyvz  sciiliwy  \v  palisadzie 
turbiny  i spr?zarki.  Badano  eksperymentalnie  i u.vzgledniono  przy  rozwirpzywaniu  rownah  rucha  cztjs'.ek 
ntechanizm  zderzertia  i odbicia  czrjstck  stalych  od  seianki  lopatki.  Badano  wplyw  sredniej  srednicy  czas- 
tek  stalych,  gesiosci  rnaterialu,  z ktorego  stj  zbudowane  oraz  poczqtkowych  predkosci  czastek  i gazu 
na  ich  charakterystyki  dynamiczne  w palisadzie.  Ponadio  pr/eprowadzono  badania  dossviadczalne  majace- 
r.a  celu  znalezienie  wplywu  k^ta  padania,  pr?dkosci  czqstek  i ich  rozmiaro'V  na  erozje  lopatek.  Przcdysku- 
towano  obserwacje  dotyczyce  erozji  lopatek. 


IIccjie;(OMnHMe  AirHaMHKti  Tnepabtx  *iacTini  n Dpo3nn  b peuieTKe 

Pe3ioMe  ■" 

ripoHopiHTca  f ac'ter  Tpaexropm!  ><  cxopoCTeS  TBepawx  «tacTnn  b cam Maesio.M  noToxe,  npOTexaiouteM 
xepe3  TypbMHHyiO  hjim  KOMnpeccopHyio  peureTxy  tipoiimneit.  SxcnepnMeHTaaibHO  iiccnertonanca  Mexa- 
iunx  yrtnpoa  x oTpa*cem<i<  Tnep^stx  xocraii  or  ctchkk  naiaTKu,  kotopmB  npunmnsctcn  no  BHHMaane 
rtp'.i  peraenatix  ypaBMenn”:  psBB/KeHna  xacTim.  HccaepioBanoCb  BfuifiHne  cpeztzero  .aaz.MeTpa  Tsepaux  xoc- 
uttt,  rmoTHOCTH  MaTepaana,  ns  xoToporo  ohh  cocionr,  a tz-kac  iiaxaaibHbt.x  sHaiemift  c so  poem  xacTiiu 
it  raja  na  nx  AaitasfHseckTie  xapaxTcpacnixH  b pemerKe.  Kposte  Toro,  nposojUfmicb  oxcnepiiMeHTajib- 
nue  itccaejoaaHaa  c tteststo  onpe^caeHna  B ih  i a h ; i a yraa  napemis,  cxopccTit  xocthu  a fix  pa3Mepoa  iia  apo- 
3;.'io  Jicna.Tox.  Ooc>K;taioica  pe3yji!,T&7bi  na6r.Kx.TCM  nil  apes  mi  rocarox. 
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II . 


CASCADE  AERODYNAMICS 


The  design  of  turbomachine  components  has  depended,  to  the 
present  day,  on  the  data  developed  during  the  design  of  similar 
previous  components.  The  designed  model  is  usually  tested, 
and  adjustments  are  made  until  the  design  conditions  are  achieved. 
Such  methods  of  trial  and  error,  when  applied  to  the  design  of 
advanced  engines,  are  costly.  The  development  of  new  theoretical 
techniques  that  could  be  directly  applied  to  achieve  the  optimum 
design  is  therefore  very  important.  Any  significant  improvement, 
on  the  present  level  of  the  engine  performance,  can  only  come 
with  a thorough  understanding  of  the  three  dimensional  flow  in 
turbomachines . 

Secondary  flow  is  a principal  phenomena  associated  with  the 
three  dimensional  flow  in  turbomachine  compressors  and  turbines. 

It  can  be  defined  as  the  difference  between  the  actual  flow, 
and  the  flow  which  would  occur  on  two  dimensional  axisymmetric 
and  meridional  stream  surfaces.  Of  the  many  factors  that 
contribute  to  the  establishment  of  secondary  flow,  end  wall 
boundary  layer  is  the  most  important.  The  interaction  of  the 
hub  and  casing  slow  moving  boundary  layer  flow  with  the  main  flow, 
which  is  turning  through  the  blades,  results  in  the  secondary 
flow'.  This  interaction  is  caused  by  the  blade  to  blade  pressure 
gradient,  the  radial  pressure  gradient,  and  the  centrifugal 
forces.  The  blade  end  clearance  and  the  relative  motion  between 
the  blade  end  and  the  annulus  walls  are  additional  factors  that 
contribute  to  the  establishment  of  the  secondary  flow.  In  the 
work  reported  here,  the  secondary  flow  resulting  from  the  inlet 
end  wall  boundary  layers  was  investigated.  The  losses  resulting 
from  the  viscous  dissipation  of  the  secondary  flow  velocities 
and  the  flow  separation  near  the  corner  between  the  end  wall 
and  the  blade  suction  surface  were  also  included  in  this  work. 

Addition  research  work  in  cascade  aerodynamics  is  reported 
in  the  following  abstracts  and  references. 
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Sheng , James,  "Sanders'  Type  Improved  Shell  Theory  for  General 
Orthogonal  Coordinates,"  Project  Themis  Report  No.  f 
October  1969.  (AD  617818) 

Abstract 

To  achieve  higher  performance  of  the  jet  engines,  the  recent 
trend  is  to  replace  the  solid  cross  section  of  the  compressor 
blades  and  turbine  buckets  by  thin  shells.  There  is  inconsistency 
in  Love's  first-approximation  theory  for  these  shells.  J.L. 
Sanders,  Jr.  made  remedies  to  these  defects  based  on  the  lines- 
of-curvature  coordinates.  However,  for  blades  and  buckets,  this 
kind  of  coordinate  system  is  not  appropriate  for  use  and  should 
be  avoided.  In  analogous  manner,  equilibrium  equations  based 
on  any  arbitrary  orthogonal  curvilinear  coordinate  system  for 
this  purpose  and  for  other  purposes  in  general  are  established. 


Wells , W.  and  Pavri , R.E.,  "Estimation  of  Supersonic  Compressor 
Efficiencv,"  Project  Themis  Report  No.  70-15,  October  1970. 

(AD  715005) 

A.bstract 

A semi-empirical  determination  of  the  adiabatic  efficiency 
of  a supersonic  compressor  stage  is  obtained.  The  method  utilizes 
the  normal  shock  relations,  a recently  developed  efficiency 
expression  for  subsonic  cascade  flow  and  experimental  pressure 
loss  data  from  supersonic  cascades  of  various  geometrical  shapes. 

A correction  factor  is  determined  which  makes  it  possible  to 
convert  the  subsonic  stage  efficiency  to  the  case  of  supersonic 
flow  involving  shock  losses. 


Hosny,  W,  and  Tabakoff,  W. , "An  Analysis  of  Losses  and  Secondary 
Plow  in  Turbine  Cascades,"  Project  Themis  Report  71-23,  December 
1971 . (AD  736853) 

Abstract 

A general  review  of  different  cascade  losses  is  presented 
with  emphasis  on  the  secondary  flow.  An  analytical  analysis 
utilizing  the  Squire  and  Winter  secondary  flow  analysis  along  with 
the  triangle  model  for  end  wall  boundary  layer  proposed  by 
Taylor  is  detailed  in  an  attempt  to  show  the  behavior  of  stream- 
lines near  the  end  wall. 
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Calculation  of  Supersonic 
Compressor  Losses 

W.  R.  Wells*  and  W.  TabakoffI 
University  of  Cincinnati,  Cincinnati,  Ohio 

Nomciiclat  urc 

.'1 1 /A ! I = S/S  t = dump  are:i  ratio 
= cliord 

= correction  factor 
= 

- radial  blade  length 
it  = Mach  number 

/’  - pressure 

S = blade  spacing 

1 = ratio  of  specific  heats 

0 = blade  camber  angle 

v --  stage  adiabatic  efficiency 

<7  - c/S  = solidity 

r - trailing  edge  thickne-s 

£ »=  cascade  stagger  angle,  measured  from  axial  direction 

Subscripts 

0 = subsonic  conditions 

1 — conditions  upstream  ahead  of  cascade  inlet 

2 ■»  condition’s  downstream  of  cascade 

n 7=  nominal  conditions 
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s = conditions  behind  shock  wave 

t — total  conditions 

I.  Introduction 

r?MlE  theoretical  prediction  of  supersonic  flow  properties  in 
-*■  a curved  channel  or  between  adjacent-  blades  of  a cascade 
presents  a formidable  task  because  of  the  complex  nature  of 
the  interaction  of  shock  waves,  the  vortex  sheets,  and  the 
boundary  layer.1-5  A more  reliable  means  of  obtaining  the 
performance  of  these  compressors  is  experimental  testing, 
as  described  in  Refs.  3-7.  For  design  purposes  the  idea  of  a 
simple  analytical  or  seinicmpirieal  means  of  predicting  the 
performance  is  attractive.  Several  papers  along  these  lines 
have  appeared  recently  in  the  literature,  among  them  is  the 
work  of  lialzer*  in  which  the  boundary-layer  blockage  effect 
and  change  in  shock  position  are  accounted  for.  A semi- 
empirical  method  for  predicting  the  performance-  of  high  re- 
aci  ion  supersonic  compressor  blade  sections  is  given  by 
Boxer.’ 

The  present,  analysis  is  an  attempt  to  supplement  the  pre- 
vious analyses  with  still  another  semiempirical  performance 
estimation  that  is  believed  to  bo,  simpler  in  application  and 
more  widely  applicable  to  a largo  family  of  cascade  geometries. 
1 ho  success  of  the  method  is  due,  in  part,  to  the  manner  in 
which  the  experimental  data  was  used  to  determine  the  initial 
Mach  number  influence  on  a key  parameter  in  the  efficiency 
expression.  This  method  extends  a successful  formulation 
developed  for  subsonic  compressors  by  Losev  and  Tabakoff*’ 
to  the  cases  of  supersonic  compressors  in  which  shock  losses 
are  present  and  accounted  for.  Other  effects,  such  as  dump 
losses  and  errors  in  the  shock  structure  model  used,  are  par- 
tially accounted  for  through  the  use  of  the  expel  imcntal  data. 

II.  Mathematical  Flow  Model 

The  purpose  of  this  analysis  is  the  development  of  a simple 
realistic  means  to  compute  the  adiabatic  efficiency  of  a supcr- 
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Weidenhamer , Gerald  H.  and  Sheng , James,  "On  the  Free  Vibration 
of  a Lone;  Helicoi dal  Membrane  Strip,"  Projec t Themis  Report  No. 
72-32,  September  1972. 

Abstract 

The  equations  of  motion  for  a shell  are  developed  in  general 
orthogonal  coordinates  using  Hamilton's  principle.  They  are 
modified  to  yield  a set  of  improved  equations  v/hich  are  then 
transformed  to  helicoidal  shell  coordinates.  The  membrane 
portion  of  these  equations  are  solved  for  the  natural  frequencies 
and  associated  mode  shapes  for  a uniformly  twisted  membrane 
boundary  valve  problem. 


Iicsnv,  W.  and  Tabakoff,  W. , "Numerical  Soluti on  for  Th ree  Dimensional 
Rotational  Flow  in  Cascades,"  Report  No.  76-48,  June  1976  . 

Abstract 

A theoretical  analysis  is  established  to  deal  with  the 
inviscid  three  dimensional  rotational  flow  which  results  from 
the  secondary  vorticity  in  turning  channels.  The  equations  of 
motion  are  written  and  manipulated  in  a special  form  to  include 
the  secondary  vorticity  as  a variable  among  the  other  flow 
variables.  The  equations  in  their  forms  are  not  attractable 
to  analytic  solutions  and  therefore  numerical  solutions  are 
adopted.  The  solution  to  the  problem  in  general  yields 
to  the  full  picture- of  the  three  dimensional  flow  together 
with  the  different  characteristics  of  the  secondary  vorticity 
distribution  in  the  flow  field.  Since  there  was  no  limitation 
on  the  flow  deflection  angle,  the  channel  radius  of  curvature 
or  the  flow  inlet  shear,  the  analysis  is  then  recommended  to 
be  used  to  obtain  the  rotational  flow  in  the  channels  of 
highly  leaded  turbomachinery  cascades. 
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JET  MIXING  FLOW 


Previous  analyses  have  shown  that  at  medium  and  low  flight 
Mach  numbers  an  exhaust  gas  jet  of  relatively  small  mass  and 
high  velocity  is  an  inefficient  method  of  producing  thrust. 

This  is  due  to  the  high  energy  losses  at  the  nozzle  exit  and 
the  poor  use  of  the  kinetic  energy  of  the  exhaust  jet  for  the 
production  of  thrust.  Two  other  important  considerations  for 
current  engines  are  good  fuel  economy  and  a low  noise  level. 

In  order  to  accomplish  these  objectives  of  maximizing  thrust, 
lowering  fuel  consumption  and  decreasing  noise,  a mixed  flow 
turbofan  engine  is  the  optimal  candidate. 

A high  airflow  is  necessary  to  increase  thrust  while  a 
high  jet  velocity  will  increase  thrust  but  decrease  propulsive 
efficiency.  The  turbofan  engine  produces  a compromise  between 
maximum  airflow  and  maximum  jet  exhaust  velocity.  Theoretical 

studies  have  shown  that  the  thrust  of  a turbofan  engine  can  be  1 

increased  by  mixing  the  hot  primary  jet  with  the  colder- 
secondary  airflow.  Many  investigators  have  confirmed  these 
theories . 

Therefore  it  is  important  to  optimize  the  amount  of  mixing 
in  a turbofan  nozzle.  In  order  to  optimize  this  mixing,  it  is 
necessary  to  better  understand  the  mechanics  of  turbulent  jet 
mixing.  Although  an  actual  turbofan  engine  may  have  turbine 
exit  swirl  and  an  augmentor  flame  holder  in  the  flow  path,  a - 
simplified  analysis,  which  is  both  reliable  and  accurate,  is 
necessary  before  the  more  complex  problem  can  be  solved.  To 
satisfy  this  need  for  a tenable  turbulent  mixing  theory,  we  have 
made  analytical  studies.  The  theory  proposed  seams  to  provide  j 

a fairly  simple,  yet  promising,  approach  to  the  turbulent 
mixing  problem.  The  theory  lacked  an  experimental  investigation 
of  the  predicted  results.  Therefore,  additional  studies  ware 

performed  to  experimentally  verify  the  accuracy  of  the  theore-  j 

tical  analysis  of  turbulent  jet  mixing  between  two  concentric 
air  streams  in  a constant  area  duct.  An  eccentric  mixing 
configuration  was  used  as  an  additional  experimental  study. 

Non-isoenergetic  turbulent  jet  mixing  experiments  were  conducted, 
where  the  measured  temperatures  were  compared  with  the  mixing 

theory.  These  two  non-isoenergetic  mixing  studies  were  believed  j 

to  have  resulted  in  a better  understanding  of  turbulent  mixing 

in  a turbofan  engine.  1 

The  research  work  performed  in  this  area  is  presented  in 
the  following  abstracts  and  references. 
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Khanna , K.K.  and  Tabakoff,  W. , "A  Study  of  Non-Isoenerqetic 
Turbulen t Jot  Mixing  Re  tween  Compressible  Subsonic  Streams  in 
Ax i- Symmetric  Constant  Area  Duct,"  Project  Themis  Report  No. 

69-1,  August  1969.  (AD-696439) 

Abstract 

The  problem  of  turbulent  mixing  between  two  compressible 
streams  in  a constant  area  duct  is  analyzed.  The  velocity  and 
temperature  profiles  have  been  worked  in  the  main  region  of  the 
mixing  chamber  for  axi~symmetric  case.  Turbulent  free  jet  experi- 
mental data  is  used,  to  obtain  expressions  for  velocity  and 
temperature  as  a function  of  the  distance  from  the  inlet  to  the 
mixing  chamber.  It  is  found  that  energy  diffuses  more  rapidly 
than  momentum  for  ducted  mixing.  The  results  of  the  theoretical 
analysis  indicate  increasing  tendency  towards  segregation  of  the 
streams  for  decreasing  initial  velocity  difference  between  the 
mixing  streams.  On  the  basis  of  the  mixing  chamber  length 
required  to  achieve  adequate  degree  of  mixing,  it  is  concluded 
that  a corrugated  rim  primary  nozzle  is  needed  to  accelerate  the 
mixing  rate  by  increasing  the  contact  surface  between  the  streams 


Ta'.alo  f _f_ , W.  _ an  <3  U o sny , _W . M.  , "The  oretica]  _a  r.d  P x peri  mental 
Jnvestiec.tions  on  the  Mixing  of  Isoenerget.ic  Confined  Co-Axial 
Jets,  ""Project  fVrpl¥~Hoport  No.  70-lQ  , ' June  1970.  Tad  710284) 

Abstract 

A.  survey  of  experimental  results  for  turbulent  mixing 
between  two  compressible  streams  in  a constant  area  duct  is 
presented.  The  experimental  results  are  for  area  ratios  3 
and  7.16,  for  velocity  ratios  ranging  from  0.4  to  2.8,  and 
with  velocity  magnitudes  varying  from  200  ft/sec.  to  950  ft/sec. 
A new  constant  of  turbulence  is  recommended  for  the  theoretical 
analysis  of  the  confined  jet  mixing  which  uses  the  free  jet 
mixing  concepts.  The  regions  in  which  this  analysis  could  be 
used  are  delineated. 


Tabakoff, W.  and  Hosnv,  W.M. , "Theoretical  and  Experimental 

Jet  Mixing  of  an  Eccentric  Primary  Jet  in a Constant  Area  Duct," 

Project  Th e m i s Report  No.  70-11,  July  1 970  . (AD  712333) 

Abstract 

The  problem  of  turbulent  mixing  between  two  compressible 
streams  in  a constant  duct  area  and  with  eccentric  primary  flow 
jet  is  analyzed.  A theoretical  analysis  is  presented  for  the 
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velocity  profiles.  It  was  found  that  these  profiles  with  or 
without  eccentricity  are  almost  the  same  if  they  are  referred 
to  the  axis  of  the  primary  flow  jet.  An  experimental  investi- 
gation for  the  case  of  eccentricity  ratio  (e/R)  equal  0.25, 
given  area  ratios  3 and  7.16  and  with  a velocity  ratio  ranging 
from  1.2  to  2.9  was  performed. 

The  theoretical  analysis  shows  good  agreement  with  the 
experimental  results,  especially  for  cases  with  high  area 
ratios  and  with  low  velocity  magnitudes. 


Ghia , K , N , , "Analytical  Investigation  of  Confined  Turbulent 
Mixing  of  Jets , " Project  Themi s Repo rt  72-24  , June  19 7 2 . 

Abstract 

The  turbulent  mixing  of  confined  axisymmetric  incompressible 
jets  was  studied,  using  the  turbulent  boundary  layer  equations  in 
terms  of  transformed  von  Mises  variables,  with  semi- empirical 
formulations  for  the  turbulent  transport  coefficients.  The 
solution  was  obtained  by  the  method  of  finite  differences. 
Analytical  predictions  of  the  jet  mixing  flow  fields  were  obtained 
for  several  configurations  and  covered  a reasonably  wide  range 
of  values  of  the  problem  parameters.  The  velocity  ratio  U^/l^ 
ranged  from  0.438  to  5.0,  the  radiu:,  ratio  Pg./R  from  0.167  to 
0.5  and  the  by-pass  ratio  X varied  from  1.626  to  16.37.  The 
average  flow  velocity  UaVq  ranged  from  abort  52  ft/sec  to 
approximately  380  ft/sec. 

Correlation  of  the  analytical  solutions  with  the  corresponding 
available  experimental  data  leads  to  the  following  principal 
results  For  confined  turbulent  jets,  the  mixing  becomes  faster 
as  the  velocity  ratio  U i/U 2 deviates  from  unity,  or  as  the  radius 
ratio  Rj/R  reduces  below  0.5,  or  as  the  average  flow  velocity 
Uavg  decreases.  The  initial  mixing  region  manifests  a rise  in 
the  local  static  pressure  and  an  accompanying  deceleration  of 
the  flow.  The  magnitude  of  this  pressure  rise  decreases  as  the 
amount  of  mixing  becomes  smaller.  A turbulent  mixing  model 
which  satisfactorily  predict;-  mean  flow  profiles  for  a reasonably 
broad  spectrum  of  flow  conditions  has  been  formulated.  The 
behavior  of  the  correlating  parameters  appearing  in  this  model 
with  respect  to  the  by-pass  ratio  X has  been  determined  and 
presented  graphically.  Analysis  of  additional  data  may  lead  to 
further  refinements  of  this  behavior. 
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Theme 

T1  IE  problem  of  turbulent  mixing  between  two  incompress- 
ible streams  in  a constant  area  duct  with  an  eccentric 
primary  How  jet  is  analyzed.  A theoretical  analysis  is  present- 
ed for  the  velocity  profiles.  An  experimental  investigation 
was  performed  for  an  eccentricity  ratio  (c  R)  of  0.25,  area 
ratios  of  3 and  7.16,  and  a velocity  ratio  ranging  front  1.2  to 
2.9. 

The  theoretical  analysis  shows  good  agreement  with  the 
experimental  results,  especially  for  the  cases  with  high  area 
ratios  and  with  low  velocity  magnitudes. 

Contents 

In  many  industrial  applications  it  is  necessary  to  deal  with 
a jet  expanding  eccentrically  into  a confined  stream  of  fluid. 
Typical  examples  arc:  jet  engines,  by-pass  ducts,  arrange- 
ments for  reduction  of  aerodynamic  noise  level  and  many 
others. 

The  theoretical  analysis  of  confined  coaxial  jet  mixing  is 
presented  by  Tabuk  iff  and  Khanna.1  In  their  analysis  they 
assumed  that  velocity  profiles  in  the  main  region  are  similar 
to  those  of  a freejet  for  which  the  universal  function  of  the 
nondimcnsicnal  excess  velocity  is 

Al/;A[/„  = (1  -r  5)1  (0 

where  £•--  Yir;  r ■ - fr-.cict  radius;  Y — distance  front  the 
axis  in  the  transverse  dhection;  AC/  = excess  velocity  at  Y\ 
and  Af/„  — excess  velocity  at  axis. 

In  the  cu  e of  coaxial  mixing,  the  velocity  profile  is  taken  to 
be  the  portion  of  the  nondinionsional  velocity  profile  which 
lies  between  Y — R and  Y - — R,  where  R is  the  radius  of  the 
duct.  For  eccentric  mixing,  the  velocity  profi'es  are  assumed 
to  be  the  eccentric  portion  of  the  freejet.  These  profiles  are 
not  symmetric  with  the  mixing  duct  axis,  however,  they  are 
symmetric  with  the  primary  How  axis.  Figure  1 shows  the 
physical  model  with  the  corresponding  notation. 

In  the  c..se  of  a freejet,  it  was  concluded  that,  owing  to  the 
similat ity  irt  velocity  profiles,  /,'r  will  remain  . onstant  (/  is  the 
mixing  length).  Moreover,  dr,\ix  — const  — C or  r = Cx. 
Now,  ~ (r)  fconst)  (Cx)(c<  int);  from  this  last  relation  it  is 
clear  that  the  mixing  length  depends  or.  C.  Since  the  mixing 
length  is  the  main  parameter  which  specifies  the  degree  of 
turbulence,  tiic  constant  C is  named  the  “turbulence  constant." 
This  constant  was  found  to  be  equal  to.  0.2  in  the  case  of  free 
jet  mixing  and  0.7  for  confined  jet  mixing.3  For  confined 
eccentric  mixing,  the  turbulence  constant  C is  taken  to  be 
equal  to  0.7. 
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■ INITIAL.  ReBCn  MAIN  REGION 


Fig.  1 Confined  jet  mixing  with  eccentricity. 


The  fundamental  equations  governing  the  flow  of  iso- 
energetic  eccentric  jet  mixing  tire  developed  as  follows.  The 
continuity  equation  for  one-dimensional  flow  for  a control 
volume  beginning  at  an  arbitrary  cross  section  at  a distance 
x,  and  ending  at  the  section  where  a uniform  flow  is  attained. 
nitty  be  written  as  follows: 

I pUdA—p3U3A3  (2) 

J O 

where  the  subscript,  3,  refers  to  the  uniform  flow  properties 
after  complete  mixing. 

In  order  to  integrate  Eq.  (2),  using  the  universal  excess 
velocity  profile  for  turbulent  freejets,  as  shown  by  Abramo- 
vich,2 one  subtracts  the  quantity 

p3U3A3~  I PU:dA 

J 0 

from  both  sides  of  Eq.  (2),  and  using  polar  coordinates  the 
following  equation  is  obtained 

A(/3p3w/?2  — I f pA  UdydO  (3) 

Jo  Jo 

where  y„  is  the  wall  radius  measured  from  the  primary  flow- 
axis. 

Dividing  each  side  of  Eq.  (3)  by  the  quantity  A U„~R2, 
introducing  f = y’r ; f,  = R;r;  {»  = yjr,  and  substituting 
Eq.  (I)  into  Eq.  (3),  for  the  incompressible  case,  the  following 
expression  is  obtained: 


Upon  substitution  of  a:  c,  (1  4-fcos 0),  where  e~c!R 
one  obtains 

MU/ML  = A,({.)  -F  (5) 

where  /?,(£, ,c)  = 0.5  — 2.5f,1-5  4-  f»5(2  -r  0.75e2)  and  .4,(1,) 
= 1 +0.4f,3-  1.143f,'-5. 
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Blasenak,  J.H.  and  Tabakoff,  W. , "Study  of  Non-isoenergetic 
Turbulent  Jet  Mixing  in  a Constant  Area  Duct,"  Report  No. 

76-49,  September  1976. 

Abstract 

A study  of  non-isoenergetic  turbulent  jet  mixing  between 
two  streams  has  been  conducted.  Using  a previously  derived 
theoretical  analysis  for  ducted  mixing,  an  experimental  investi- 
gation was  performed  to  verify  this  theory  and  to  determine 
the  non-isoenergetic  turbulent  jet  mixing  characteristics  in 
a constant  area  duct.  Temperature  profiles  were  measured  at 
several  axial  locations  in  the  duct  for  both  a concentric  and 
an  eccentric  configuration.  It  w as  determined  that  the 
theoretical  and  experimental  temperature  profiles  agreed  fairly 
well  for  both  cases,  although  the  concentric  case  showed  better 
agreement  than  the  eccentric  case.  It  was  also  determined  that 
a new  constant  of  turbulence  in  the  initial  region  was  needed 
for  non-isoenerget i c mixing,  mixing  is  generally  more  rapid  than 
the  theory  predicted,  the  initial  temperature  difference  between 
the  two  streams  did  not  have  much  effect  on  the  rate  of  mixing 
and  a higher  area  ratio  produced  better  agreement  between  the 
theory  and  the  experimental  data.  It  was  concluded  that  the 
theory  was  good  for  a fairly  simplified  analysis. 
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IV. 


HIGH  TEMPERATURE  MATERIALS 


The  high  temperature  requirements  of  future  advanced  pro- 
pulsion systems  necessitates  the  careful  consideration  of 
materials  to  be  used  in  turbine  blades,  nozzles,  and  other 
high  temperature  components.  The  absolute  upper  limit  of 
operating  temperatures  of  metal  components  of  an  engine  is 
dictated  by  the  solidus  temperatures  of  the  alloys  used. 

Through  the  phenomenon  of  dispersion  hardening  the  design 
engineer  can  develop  alloys  that  have  good  strength  and  creep 
resistance  at  appreciable  fractions  of  the  solidus  temperature. 
In  addition  to  strength  at  higher  temperatures,  suitable  jet 
engine  alloys  must  be  resistant  to  attack  by  oxidation  and 
corrosion  from  the  engine  gases.  Superalloys  based  on  nickel 
with  major  amounts  of  chromium  and  cobalt,  and  lesser  amounts 
of  titanium,  aluminum,  and  molybdenum  give  workable  alloys  of 
good  strength  at  high  temperatures,  and  good  corrosion  and 
oxidation  resistance.  Through  constant  efforts  the  operating 
temperature  of  alloys  has  been  slowly  increased  over  the  years. 
Columbium  metal  is  a potential  high-temperature  construction 
material  for  gas-turbine  buckets,  because  it  retains  a useful 
degree  of  strength  at  higher  temperatures  where  conventional 
creep  resisting  alloys  prove  to  be  inadequate.  The  major 
obstacle  in  using  columbium  in  pure  or  unprotected  form  at  high 
temperatures  is  its  poor  oxidation  resistance  in  oxygen-bearing 
environments.  Alloying  with  other  elements  has  already  proved 
a partial  rectification  of  this  disadvantage  and  although 
considerable  knowledge  about  the  oxidation  behavior  of  these 
alloys  has  been  accumulated,  both  from  basic  and  engineering 
studies,  reaction  mechanism(s)  are  not  clearly  understood. 
Principal  hinderance  in  interpreting  the  alloy  oxidation  is  the 
insufficient  information  on  the  oxidation  of  pure  columbium 
because  of  complexity  of  the  oxidation  rate-temperature 
relationship. 

The  main  purpose  of  this  study  was  to  investigate  in  detail 
the  oxidation  behavior  of  Cb-10  at  . %W  alloy.  To  supplement  the 
understanding  of  the  oxidation  mode  of  Cb-10  at  . %W  alloy, 
attempt  has  also  been  made  to  resolve  the  unusual  oxidation 
rate-temperature  dependence  of  the  columbium  metal.  The  experi- 
mental work  was  essentially  confined  to  oxidation  over  a 
temperature  range  of  900-1200 °C  in  dry  oxygen  at  100-760  torr. 


The  research  work  performed  in  this  area  is  reported  in 
the  following  abstracts  and  references. 


Sikka , V.K.  and  Rosa,  C.J.,  "Determination  of  Oxygen  Diffusion 
Coefficients  in  Tungsten  Oxide,"  Project  Themis  Report  Mo. 
70-13,  August  1970.  " (AD  7 1 500 6 ) 

Abstract 

"Interruption  Kinetic  Technique"  experiments  carried  out 
on  tungsten  oxide  films  formed  on  tungsten  yield  the  following 
experssion  for  the  diffusion  coefficient  of  oxygen  anion 
vacancies : 

D = 6.83  x 10~2  exp  (-29 , 889/RT) 
for  the  temperature  range  568°  - 908°C. 

Calculations  are  made  for  the  concentration  of  oxygen 
vacancies,  for  their  free  energy  of  formation  and  for  ionic 
conductivity  in  tungsten  oxide.  Deviations  from  stoichiometry 
in  the  oxide  are  plotted  against  temperature  on  the  available 
W-0  phase  diagram. 


Sikka,  V.K.  and  Rosa,  C.J.,  "High  Temperature  Oxidation  of 
Columbium  and  CB-10  At  . %W  Alloy,"  Project  Themis  Report  No. 
70-14,  September  3 970  . (AD-  715r007l 

Abstract 

Oxidation  kinetics  of  Cb  and  Cb-10  at  . %K  alloy  over  the 
temperature  900-1200°C,  and  oxygen  pressures  of  100-760  torr 
were  determined  using  highly  sensitive  Cahn  Electrobalance. 

The  products  of  oxidation  were  examined  metallographically 
and  by  x-ray  diffraction  methods. 

Activation  energy  for  the  parabolic  rate  constant  has 
confirmed  a diffusion  controlled  growth  process  for  the  oxidation 
of  Cb-10  at  . %W  alloy.  The  overall  oxidation  kinetics  of  Cb-10 
at  . % alloy  in  the  temperature  range  1100-1200°C  is  best 
described  by  the  paralinear  model. 

Possible  explanations  for  the  abnormal  oxidation  rate- 
temperature  relationship  for  the  oxidation  of  Cb  and  Cb-10 
at  . t,w  alloy  have  been  discussed. 
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A theoretical  solution  by  finite  difference  methods  is  presented  for  the  steady  state  temperature  .d 
with'n  a solid  uniform  ri^ht  circular  cylinder  which  is  healed  by  a unin. mi  souree  assn;  it-  ear  ed  surface 
and  which  simultaneously  loses  heat  by  fourth-power  radiation  front  all  its  senates-  to  touch  coolei  sur- 
rounding - 'Ihe  results  of  the  solution  are  more  accurate  than  the  tcmi-enipi: leal  relari  of  flush 

tt  0 1.,  and  so  tan  be  applied  with  greater  reliability  to  their  new  metal  envelope  induction  tr  r.nsp.e  for 
determining  thermal  conductivities  of  nounietals  at  high  temperatures. 


INTRODUCTION 

Several  years  ago,  I loch  el  al.'  introduced  an  experi- 
mentally simple  method  for  indirect iy  measuring  the 
thermal  conductivity  of  certain  it.  iterials  at  high  tem- 
perature. In  this  method,  a solid  right  circular  cylinder 
of  the  material  is  placed  coaxially  within  an  induction 
coil  and  heated  electromagnetically  in  vacuum.  At 
steady  state,  the  radial  tempera' lire  profile  a;  one  o: 
the  flat  faces  of  the  cylindei  is  me  isured  with  an  optical 
pyrometer.  Then,  combining  these  measurements  with 
the  emissivity  and  the  dimensions  of  the  cylinder,  the 
thermal  conductivity  is  calculated. 

'-  ,e  method  has  been  used  for  Mo,1-  Va,1  and  graph- 
ite.’ Recently,  the  method  was  extended  to  materials 
such  as  Al=Oj  by  covering  the  curved  surface  of  the 
cylinder  with  a metal  envelope.'- 

Tor  cither  technique,  the  calculation  of  /.-  i-  based  on 
a highly  simplified  tin  :ma!  model.1  It  assumes  that  the 
curved  surface  is  isothermal,  and  that  the  temperature 
at  the  flat  surfaces  varies  quadraticallv  with  radial  dis- 
tance. These  assumptions  have  been  shown  to  be  in- 
accurate when  applied  to  the  original  technique  of  direct 
induction.5  The  demonstration  was  bv  means  of  a 
theoretical  solution  of  the  electric  field,  followed  by  a 
theoretical  solution  for  the  thermal  field.  However, 
based  on  these  more  accurate  results,  a simple  and  more 
reliable  method  for  calculating  k was  presented.5  * 

The  objective  of  the  present  paper  is  to  examine  the 
validity  of  Hoch’s  assumptions  as  applied  to  the  later 
envelope  induction  technique  and  to  improve  on  thei  . 
as  necessary.  This  examination  will  be  conducted  by 
employing  more  realistic  assumptions,  and  then  solving 
for  the  temperature  field  numerically  via  finite  differ- 


ij"  rate  of  heat  input,  per  unit  surface 
r radial  distance 
t absolute  temperature 

absolute  temperature  ah  :;g  the  circular  edge  <r  a, 

•=L) 

le  absolute  temperature  at  the  center  of  the  flat  fact 
(r  = 0,  z—  L) 

t„  absolute  temperature  air  eg  the  center  of  the  circa!.  : 
face  (r=  c,  z-—Q) 
z axial  distance 

a total  emissivity 

t Stefan- Iloltzmann  radiation  constant 
THEORY 

The  cylinder  is  shown  in  Fig.  1.  By  virtue  of  angular 
symmetry,  Fourier’s  conduction  equation  for  a material 
of  isotropic  constant  k at  steady  state  with  no  io!un;elric 
heat  sources  can  be  written  as  follows: 

(0-1 /Or)  + r~l  (01/ dr)  + (o-l/dz-)  = 0.  ( 1 ) 

Equation  (1)  is  the  governing  differential  equation. 

Symmetry  also  gives  Eqs.  (2),  (3)  as  two  of  the 
boundary  conditions. 


NOMENCLATURE 


a radius  of  cylinder 
k thermal  conductivity 
L half-length  of  cylinder 


tW(0,s)/dr=0  (2) 

0t(r,  0)/<te~0.  (3) 

For  convenience,  attention  can  now  be  restricted  to  the 
first  quadiant  in  Fig.  1. 

The  cylinder  is  assumed  to  be  much  hotter  than  its 
surroundings,  so  that  thermal  radiation  from  the  sur- 
roundings to  the  cylinder  is  negligible  compared  to 
that  from  the  cylinder.  A conduction-radiation  balance 
at  a differential  area  element  in  the  flat  surface  yields 
Eq.  (4)  as  the  third  boundary  condition. 

dl(r,  L)/dz=  - (icr/k)t*(r,  L)  (4) 

Due  to  the  vacuum,  there  is  assumed  to  be  no  conduc- 

tion between  the  cylinder  and  its  surroundings. 

The  heat  which  is  furnish  'd  by  the  metal  envelope 
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Abstract 

Oxidation  kinetics  of  Cb-10  at  . % Cr  alloy  were  studied  at 
three  temperature  levels:  900,  1000  and  1100°C  under  oxygen 

pressures  ranging  from  100-760  torr.  At  1000  and  1100°C  the 
oxidation  kinetics  are  lower  than  for  pure  columbium.  This 
improvement  of  oxidation  resistance  of  the  alloy  is  mainly 
attributed  to  the  modifying  effects  of  Cr  cations  in  a-Cb^-O,. 
by  stabilizing  the  compact,  inner,  oxide  layer.  The 
reaction  mechanism  is  essentially  inward  oxygen  diffusion 
across  the  compact  inner  scale.  At  900 °C  the  oxidation  kinetics 
of  the  alloy  are  higher  when  compared  with  pure  columbium  and 
follow  a linear  rate.  This  decrease  in  oxidation  resistance  is 
related  to  polymorphic  transformation  in  Cb20g  which  leads  to 
fracture  of  the  highly  distorted  a-Cb2~05.  The:  reaction 
mechanism  at  S00°C  could  be  associated  with  oxygen  absorption. 
X-ray  diffraction  patterns,  electron  microprobe  analyses  and 
the  observed  morphologies  of  the  oxide  scale  support  the 
proposed  oxidation  mechanisms  of  the  alloy. 


Victh , D.  L.  and  Pool  , M . J.  , "Desi on  and  Con s_t rncti on  of  a H i ch 

Temperature Liquid  • t -j  Coluticn  Calorimeter , " Project  ‘iixr.is 

Report  No.  72-30,  October  1972 . 

Abstract 

The  design  consideration  and  calibration  of  a new  single 
well  liquid  metal  solution  calorimeter  for  use  at  high  temperature 
is  described.  The  calorimeter  has  an  operational  temperature 
range  of  500°C  to  1300°C  in  a vacuum  and  was  designed  primarily 
to  operate  1200°C  using  liquid  uranium  as  a solvent. 

Calibration  of  the  instrument  indicates  that  the  energy  equivalent 
ranged  from  0.8100  joules/yV  (.1950  cal/pV)  to  0.6950  joules/pV 
(0.160  cal/uV) . Calibration  was  based  on  the  heat  content, 
between  room  temperature  and  the  calorimeter  temperature,  of 
weighed  samples  of  pure  uranium.  The  partial  heat  of  solution 
of  vanadium,  chromium,  iron,  and  nickel  in  uranium  at  1150 °C 
were  determined. 
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Birla,  N.C.  and  Hoch,  M. , "The  Solutioning  and  Aging  Reactions 
and  Thoi  r h'f  f< -ct  < *n_  Hardness  and  M icrost  rue  tun  o in  Borne 
Co  1 u.iib  i r 1 m Base  A1 1 oys,"  ?-*roject  Themis  Report  No.  72-31, 

September  1 9 7_2  . 

Abstract 

A study  was  made  to  investigate  the  solutioning  and  aging 
reactions  and  their  effects  on  microstructure  and  hardness  in 
twenty-four  experimental  columbium-base  alloys,  involving  Cb-C, 
Cb-Si-C,  Cb-15Hf-C,  Cb-15Hf-C-Si , Cb-33IIf-C  and  Cb-33Hf-C-Si 
systems.  After  selected  heat  treatments,  the  microstructural 
changes  were  determined  by  electron  microscopy  and  the  phases 
were  extracted  and  identified  by  x-ray  diffraction.  This  study 
involved  more  than  72  extractions  and  their  analysis  by  x-ray 
diffraction,  more  than  135  solution  treatments  in  an  induction 
furnace,  about  75  heat  treatments  for  aging  in  a Brew  furnace 
and  more  than  550  samples  for  polishing,  metallography  and 
microhardness  measurements . 

The  aging  reactions  in  the  Cb-15Hf-C,  Cb-15Hf-C-Si , and 
Cb-33Hf-C-Si  alloys  are  characterized  by  precipitation  of  the 
complex  monocarbide  (Cb,  Hf)  ( C , O,  N)  phase. 

Replacement  of  some  of  the  carbon  by  silicon  in  the  Cb-15Hi-lC 
alloys  raises  the  overaging  temperature  from  80CQC  to  1000°C  and 
leads  to  strength  retention  for  longer  times  than  the  alloys  which 
contain  carbon  only.  However,  the  alloys  in  Cb-33Hf-C  system, 
did  not  show  precipitation  hardening  due  to  the  precipitation  of 
coarse  Cb2C  particles  during  aging  process. 
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V * GAS  TURBINE  BIADKS  AND  COMBUST  ION  CH  AM]  KPS  HKAT  TPA^NSFER 

The  efficiency  of  gas-turbine  engines  is  highly  dependent  on 
the  turbine  inlet  temperature  of  the  working  fluid.  If  the 
present  achievable  inlet  gas  temperature  is  doubled,  this  means 
theoretically,  that  a turbine  can  double  its  power  per  pound  of 
working  fluid.  As  a consequence,  the  turbine's  bulk  and  weight 
can  be  reduced.  Because  of  today's  metallurgical  limitations 
of  temperature  with  respect  to  the  blade  load  force,  production 
engines  are  limited  to  2200°F.  With  this  temperature  limitation 
on  the  blades,  the  only  way  to  increase  the  turbine  inlet  gas 
temperature  is  to  use  some  method  to  cool  the  turbine  blades. 

The  requirement  for  this  cooling  is  likely  to  remain  for  a long 
time  to  come,  even  with  the  new  materials  now  being  developed 
for  turbine  blade  use. 


The  reported  basic  research  program  was  to  investigate  the 
heat  transfer  phenomena  associated  with  internally  cooled  blades, 
namely  by  jets  impinging  and  transpiration.  A jet  impinging 
on  e surface  an  be  utilized  to  achieve  high  performance  heat 
transfer  configurations.  The  surface  may  be  cooled  by  the  jet 
flow  to  obtain  high  local  heat-transfer  directly  under  the  jet, 
while  good  overall  heat  transfer  can  also  be  obtained  when  the 
surface  area  to  be  cooled  is  appreciably  larger  than  the  jet. 

The  transpiration  cooling  uses  cooler  fluid  which  is  injected 
through  the  porous  wall  into  the  boundary  layer  of  the  blades 
to  protect  the  external  hot  surfaces.  Such  type  of  cooling  is 
well  known  to  have  high  cooling  efficiency  and  excellent 
shielding  from  severe  thermal  environment.  Critical  temperatures 
and  thermal  environments  are  encountered  in  combustion  chambers 
and  over  gas  turbine  blades.  Transpiration  cooling  perhaps 
finds  its  most  important  applications  in  keeping  these  components 
at  an  acceptable  temperature.  To  insure  efficient  performance 
over  a reasonable  lifetime  makes  the  problem  of  blade  and 
combustion  chambers  cooling  absolutely  necessary. 

The  research  work  performed  in  the  above  mentioned  areas 
is  reported  in  the  following  abstracts  and  references. 


Clevenger,  W.  and  Tabakoff,  V.’.,  "Investigation  of  the  Heat 
Transfer  Characteristic.-,  of  a Two- Dimensional  Jet  1 mninging 
on  a Semi-Cvlinder , " Project  Themis  Report  No.  69-5,  October 
1969 . (AD  697165) 

Abstract 

Boundary  layer  methods  are  ufsed  to  find  the  heat  transfer 
rate  caused  by  a two-dimensional  jet  impinging  on  the  inside 
surface  of  a semi-cylinder.  Because  of  a wider  pressure 
distribution  on  the  semi-cylinder,  the  solution  predicts  a 
delayed  transition  from  a laminar  to  a turbulent  boundary  layer. 
This  delayed  transition  causes  the  average  heat  transfer  from 
the  semi-cylinder  to  be  less  than  the  average  heat  transfer 
from  the  flat  plate.  Experimental  data  supports  the 
theoretically  predicted  heat  transfer  rate  from  the  stagnation 
area  of  the  semi-cylindrical  plate  and  indicates  heat  transfer 
rates  that  are  less  than  those  predicted  by  theory  in  the  other 
regions  of  the  flow  field. 


Tsuei,  Y . G ■ and  A3- 1 is__, E . W , , "Laminar  Jet  Mixing  with  Tnit j al 

and  Boundarv  Effects , " Project  Themis  Report  No.  70-9,  April 
1970.  (AD  704983) 

Abstract 

The  objectives  of  this  report  are  to  investigate  the  flow 
field  and  the  wall  shear  stress  for  laminar  tangential  injection. 
The  effects  of  the  boundary  and  initial  velocity  profiles  under 
zero  pressure  gradient  are  considered.  A linearized  approxima- 
tion with  a coefficient  chosen  from  the  parameters  of  the 
boundary  layer  has  been  used  to  analyze  the  slot  injection.  It 
is  found  that  wall  shear  stress  is  almost  independent  of  the 
ratio  of  jet  velocity  to  free  stream  velocity  and  the  ratio  of 
slot  height  to  initial  outside  boundary  layer  thickness  up  to 
a downstream  distance  equal  to  0.036  ReL  times  the  slot  height. 

It  is  also  found  that  the  effectiveness  of  the  jet  mixing  is 
greater  for  large  ratios  of  slot  height  to  initial  outside 
boundary  layer  thickness  than  for  small  ratios.  The  asymptotic 
solution  of  the  velocity  distribution  far  downstream  is  compared 
with  the  Blasius  profile  and  good  agreement  is  found 

The  results  can  be  used  for  film  ccoling  to  estimate  the 
spacing  between  the  injection  slots  for  effective  protection. 
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R avuri , R.  and  Tabakoff,  W. , "Heat  Transfer  Characteristics  of 
a Row  o f Air  J ots  Impinging  on  the  Jnr ide  of  a Sem-Ci  rcular 
Cylinder,"  Project  Themis  Report  No.  71-18.  January  1971. 

Tad  717111) 

Abstract 

Experimental  results  of  average  heat  transfer  coefficients 
for  a row  of  air  jets  impinging  on  the  inside  of  a semi-cylinder 
are  presented.  The  variation  of  the  average  Nusselt  Number  with 
the  wide  range  of  Reynolds  Numbers  are  found.  Two  semi-cylinder 
concave  surfaces  with  different  diameters,  four  different 
impinging  jet  configurations  and  variable  impinging  heights 
were  investigated.  The  effect  of  introducing  solid  particles 
into  the  cooling  air  on  the  heat  transfer  performance  was 
studied . 


Tabakof  f , W.  , Pavri  R . and  Cl  c venger,  W.  , "lie at  Transfer  by  a 

Multiple  Array  of  Pound  Jets Impinging  Perpendicular  to  a 

Concave  Surface,"  Project  Tnemis  Report  No.  71-20,  June  1971. 
(AD  728489)  ' 

Abstract 

An  experimental  investigation  was  undertaken  to  study  jet 
impingement  cooling  of  a semi-cylindrical  concave  surface  in  a 
semi-enclosed  environment.  The  investigation  showed  that  it 
was  possible  to  obtain  a correlating  formula  for  the  various 
parameters  involved.  The  results  are  compared  with  tht>se  of 
Kercher  and  Tabakoff  who  carried  out  a similar  investigation 
using  a flat  plate. 


Pavri,  R.  and  Tabakoff,  W,  , "An  Analytical  Solution  of  Wall 

Temperature  D i st ri but!  o n f qr_  Tr ansp i^ati on  and  Local  Mass 

Injection  Over  a Flat  Mate,"  Project  Tnemis  Report  No.  72-25, 
March  1972. 

Abstract 

This  analysis  describes  an  analytical  solution  of  the  non- 
similar laminar  boundary  layer  with  continuous  or  local  injection 
and  variable  wall  temperature.  The  solution  is  for  a flat  plate 
with  zero  pressure  gradient,  however,  it  can  be  modified  for 
variable  pressure.  The  method  consists  of  transforming  the 
partial  differential  equations  for  momentum  and  enthalpy  and 
then  solving  the  transformed  equations  by  assuming  polynomial 
velocity  and  temperature  profiles. 

The  analysis  is  presented  in  two  parts.  The  first  part 
presents  the  solution  for  continuous  injection  (transpiration 
cooling)  while  the  second  part  involves  discontinuous  boundary 
conditions  due  to  local  injection. 
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Pavri,  R.E.  an  d Tabakoff,  W.  , "An  Analytics  1 So  ] ut.lon  of  the 
Compressible  Laminar  Momentum  and  Thermal  Boundary  Layers  wl th 

P r essure  Gradient-  and  Con t.inuous  Mass  Injection L" P r o jc c h 

Themis  Report  No.  72-29,  September  1972.  (AD  751804) 

Abstract 

The  study  describes  an  analytical  solution  of  the  nonsimilar 
laminar  boundary  layer  with  pressure  gradient,  variable  wall 
temperature  and  continuous  injection.  The  method  consists  of 
transforming  the  partial  differential  equations  for  momentum 
and  enthalpy  and  then  solving  the  transformed  equations  by 
assuming  polynomial  stream  function  and  enthalpy  profiles. 
Solutions  obtained  show  very  good  agreement  with  exact  numerical 
results . 

The  solutions  are  obtained  for  flows  over  wedges  as  well 
as  at  the  two-dimensional  stagnation  point  and  over  curved 
surfaces  of  a two-dimensional  body  in  crossflow.  The  results 
of  the  study  show  that  the  boundary  layer  is  very  strongly 
effected  by  the  injection  mass  flow  rate. 
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| Gas  Turbine  Blade  Heat  Transfer 
| Augmentation  by  Impingement  of  Air  jets 
''  Having  Various  Configurations' 


An  experimental  investigation  of  heal  transfer  characteristics  for  various  configurations 
of  air  jets  impinging  on  the  leading  edge  inner  surface  of  the  blade  wall  is  presented. 
Three  configurations  'acre  investigated,  namely  a slot  jet,  a round  jet  row  and  an  array 
of  round  jets.  The  effect  on  the  heat  transfer  coefficient  of  injecting  solid  particles  into 
the  air  flow  is  considered..  The  study  treats  an  important  class  of  turbine  blade  cooling 
for  which  small  cooling  mass  f.ow  rales  arc  of  interest.  The  experimental  facility  and 
procedure  arc  described  in  detail.  A theoretical  t- : hi.tque  is  introduced  for  predicting 
the  hen!  transfer  in  the  case  of  the  clot  jet  configuration.  The  results  ore  compared  to 
experimental  data. 


1 


Introduction 

The  efficiency  of  gas-turbine  engines  is  highly  de- 
pendent  on  tire  turbine  inlet  temperature  of  the  working  fkiid. 
In  today’s  operating  gas  turbine  the  inlet  gas  temperature  is  be- 
low the  adiabatic  flame  temperature  of  hydrocarbon  fuels.  If 


1 This  work  was  sponsored  under  Project  Themis  Contract  Number 
DAHC04-G0C-00IG,  Arm;.-  Research  Office,  Durham. 

Contributed  by  the  Gas  Turbine  Division  and  presented  at  the 
Gas  Turbine  Conference  and  Products  Show,  Houston,  '1  ex.v,  March 
28-April  1.  1971,  of  The  American  Society  of  Mechanical  Engi- 
neers. Manuscript  re  eived  at  AS.ME  Headquarters,  December  21, 
1970.  Paper  No.  71-GT-9. 

Nomenclature 

Ap  - surface  area  of  the  heater 
plates,  in.5 

FS  = constant  used  to  correctly 
position  the  velocity  decay  of 
a turbulent  wall  jet,  in. 

Co  *»  constant  used  in  the  velocity 
decay  equation  of  a turbu- 
lent wall  jet 

Cp  = specific  heat  at  constant  pres- 
sure, Btu/lbm  deg  U 

I)  - round  jet  diameter,  in. 

II  *=  height  of  the  nozzle  exit  above 
the  plate,  in. 

h = local  heat  transfer  coefficient, 

Btu/hr  ft*  deg  F 

h — average  height  transfer  coef- 
ficient, Btu/hr  ft5  deg  F 

k = thermal  conductivity  of  air, 

Btu/hr  ft  deg  F 


the  present  inlet  gas  temperature  of  2,000  deg  F is  doubled,  the- 
oretically the  turbine  can  double  its  power  per  pound  of  working 
fluh].  As  a consequence,  the  tin* bi;./s  bulk  ut  1 weight  can  be 
reduced.  Because  of  today’s  metallurgical  limitations  of  tem- 
perature with  respect  to  the  blade  load  force,  production  engines 
are  limited  to  2,300  deg  F.  With  this  temperature  limitation  on 
the  blades,  the  only  way  to  increase  the  turbine  inlet  gas  tem- 
perature is  to  use  some  method  to  cool  the  turbine  blade.  In 
order  to  design  and  develop  an  effective  cooling  means,  a know- 
ledge of  the  blade  beat  transfer  i-  essential. 

The  purpose  of  this  project  was  to  compare  the  effectiveness  of 
three  different  systems  of  air  jets  impinging  on  ti  e inside  surface 
of  a half-circular  cylinder.  Such  a configuration  would  be  similar 
to  that,  of  a turbine  blade  cook'd  by  gas  jets  impinging  on  its  in- 


L = length  of  the  semicylinder 
model  and  the  length  of  the 
slot  jets,  in. 

m = exponent  used  with  the  Fraud tl 
Number  m = */* 

rii  — total  air  mass  flow  from  all  jets, 
Ibm/hr 
n±_  1 
2 

exponent  used  in  the  velocity 
and  tempemture  profiles  of 
the  wall  jet 

Nusselt  number  based  on  plate 
length  Nu  = FiL/k 
Nu.j  — local  Nusselt  Number  based  on 
the  jet  half  width  at  half 
depth  Nu.s  = hf/k 

F(x*)  « pressure  at  any  point  x*  along 
the  plate,  psia 


A'  = 


Nu  - 


P(o)  — pressure  at  the  point  x*  = 0, 
the  stagnation  pressure;  psi:i 
P»n> t = ambient  pressure,  psia 
Fr  = Frandtl  number  Fr  *=  f.iCp/k 
Qo  — gross  power  input  to  the  plate, 
Btu/hr 

Ql  = power  losses  through  the  insu- 
lated box  supporting  the 
model,  Btu/hr 

<7u/  — theoretical  heat  transfer  rate 
from  the  plate,  Btu/hr  ft5 
Re  *=  Reynolds  number  Re  = rii/pL 
St  = Stanton  number  St  ~ h/CpUjp 
s - spacing  between  tl.c  center 
line  of  the  round  jets,  in. 

Tw  = plate  wall  temperature,  deg  F 

(Continued  on  next  page ) 
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This  analysis  describes  an  analytical  solution  of  the  nonsimilar  laminar  boundary 
layer  with  continuous  or  local  injection  and  variable  wall  temperature.  The  solution 
is  for  a flat  plate  with  zero  pressure  gradient,  he  vever,  it  can  be  modified  for  vari- 
able pressure.  The  method  consists  of  transforming  the  partial  differential  equa- 
tions for  momentum  and  enthalpy  and  then  solving  the  transformed  equations  by 
assuming  polynomial  velocity  and  temperature  profiles.  The  analysis  is  presented 
in  two  parts.  The  first  part  presents  the  solution  for  continuous  injection  (transpira- 
tion cooling;)  the  second  part  involves  discontinuous  boundary  conditions  due  to 
local  injection. 
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OPTIMUM  HEAT  TRANSFER  CHARACTERISTICS  OF  SEMI-CIRCULAR 
SURFACES  COOLED  RY  AIR  IMPINGEMENT  FROM  AIR JET 
ARRAYS  AND  ROW  OF  AIR  JET  NOZZLES* 

by 

W.  TABAKOFF  and  R.  STHENECK 
Department  of  Aerospace  Engineering 
University  of  Cincinnati 
Cincinnati,  Ohio  U.S.A. 


The  trend  in  today's  gas  turbine  engine  design  is  toward 
higher  turbine  inlet  temperature.  As  a consequence  the  turbine's 
bulk  and  weight  can  be  reduced.  Because  of  today's  metalur gical 
limitation  of  temperature  with  respect  to  the  blade  load  forces, 
the  only  way  to  increase  the  inlet  gas  temperature  is  to  use  some 
method  to  cool  the  turbine  blade.  The  authors  investigated  dif- 
ferent systems  of  impinging  air  jots  cooling  the  inside  surface 
of  a -half-circular  cylinder.  Such  a configuration  will  be  similar 
to  that  of  a turbine  blade  cooled  by  gas  jets  impinging  on  the 
leading  edge  inside  surface.  This  paper  is  a continuation  of 
the  above  mentioned  work,  the  purpose  of  which  is  to  derive  general 
equations  or  scaling  factors  for  optimum  cooling  of  semi-cylinder 
surfaces  of  any  diameter  by  air  impingement  from  air  jet  arrays 
and  nozsle  rows.  In  order  to  obtain  experimental  data  5 inch, 

2.5  inch,  1.25  inch  and  0.5  inch  diameter  semi-cylindrical  models 
were  investigated. 
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Ravuri , R.  and  Tabako f f , W .,  "A  Numeric al  Solution  for  the  Heait 
Transfer  nctwcon  an  Axi -Symmetric  Air  Jot  and  a Hunted  Plate,"  ” 
Report  No.  V 3-38 , September  1973 . 

Abstract 

Due  to  their  characteristic  high  rate  of  heat  transfer, 
impinging  air  jets  are  being  used  extensively  in  industry.  A 
numerical  scheme  is  devised  to  predict  the  heat  transfer  betv/een 
an  axi-synmetric  air  jet  and  a heated  plate,  neglecting  the 
effect  of  curvature,  compressibility,  and  turbulence.  The 
momentum  and  continuity  equations  are  transformed  to  vorticity 
and  stream  function  equations  and  are  then  solved  by  iterative 
successive  substitution  techniques  to  determine  the  flow  field 
and  subsequently  the  temperature  distribution  in  the  fluid. 

An  attempt  is  made  to  predict  the  average  heat  transfer 
coefficient  of  a row  of  axi-symmetric  air  jets  impinging  on  a 
heated  plate  by  assuming  that  each  jet  cools  a fraction  of  the 
plate  without  any  interference  from  the  neighboring  jets. 

\ 


Nilson,  R.li.  and  Tsuei,  Y.G.,  "A  Numerical  Method  for  Boundary 
buyer  Donations, " Report  No.  74-41,  March  1974. 

Abstract 

A general  marching  procedure  for  numerical,  solution  of'  the 
two-dimensional  boundary  layer  equations  is  presented.  Although 
the  present  method  traces  its  origin  to  that  of  Patankar  and 
Spalding,  major  modif ications  have  been  undertaken.  The  rrethod 
is  tested  for  a variety  of  laminar  flow  configurations  including 
wedge  flows,  incompressible  flow  on  a flat  plate  at  various 
Eckert  numbers,  compressible  flow  on  an  adiabatic  flat  plate 
at  various  Mach  numbers,  flow  over  a cylinder,  uniform  suction 
on  a flat  plate,  and  Howarth  flow.  The  results  show  excellent 
agreement  with  solutions  by  other  methods,  even  when  as  few  as 
fifteen  grid  lines  are  used.  A separate  report  entitled 
"Film  Cooling  by  Oblique  Slot  Injection"  presents  results  for 
a variety  of  injection  configurations  which  show  the  effects  of 
coolant  mass  flow  rate,  injection  angle,  boundary  layer  thickness, 
slot  width,  and  multiple  slots. 
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Film  Cooling  by  Oblique  Slot  Injection 

R.  H.  Nilson*  and  Y.  G.  Tsuei+ 

University  of  Cincinnati,  Cincinnati,  Ohio 

Introduction 

\"57ALL  cooling  by  boundary-layer  injection  or  transpiration 
y V is  used  in  engineering  applications.1  The  present  in- 
vestigation concerns  the  film  cooling  effectiveness  of  oblique 
injection  from  the  wall  into  a compressible  laminar  boundary 
layer  through  single  or  multiple  slots.  Numerical  solutions  of 
the  boundary-layer  equations  are  obtained  by  a finite-difference 
method  which  has  been  extensively  tested  and  found  to  be 
accurate,  versatile,  and  very  stable.  Film  cooling  effectiveness  is 
presented  for  a wide  variety  of  injection  configurations  so  that 
the  effects  of  coolant  mass  flow,  injection  angle,  boundary-layer 
thickness,  slot  width,  and  the  presence  of  upstream  cooling  slots 
can  be  investigated.  The  results  arc  interesting,  and  conclusions 
heretofore  unreported  are  drawn  regarding  selection  of  film 
cooling  parameters. 


Received  November  6,  1973;  revision  received  December  26.  1973. 
This  work  was  partially  supported  by  U.S.  Army  Research  Office  — 
Durham,  under  Contract  DAHC  0-1-69-C-0016. 
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Tsuei , Y . G . and  Nilson,  R.H.,  "Influence  of  Injc-  cti  on  Angle  on 
Film  Cool  in g Efiec tivene ss_  i n Laminar  Compre ssible  Flow, " 

Feport  No.  75-45,  May  1975 . 

Abstract 

* 

Wall  cooling  effectiveness  is  investigated  for  tangential, 
inclined,  and  normal  injection  of  coolant  through  single  or 
multiple  wall  slots  into  a laminar  compressible  boundary  layer. 
Numerical  solutions  of  the  boundary  layer  equations  are  obtained 
by  a finite  difference  method  which  has  been  extensively  tested 
and  found  to  be  accurate,  versatile,  and  stable.  A grid  control 
procedure  which  maintains  a constant  flow  rate  between  grid 
lines  is  found  to  be  well  suited  to  the  present  injection 
calculations  wherein  the  boundary  layer  growth  in  the  slot  is 
as  much  as  hundred-fold  and  the  longitudinal  component  of  the 
injection  velocity  is  in  some  cases  as  large  as  the  free  stream 
velocity.  Film  cooling  effectiveness  is  reported  for  a variety 
of  injection  configurations  so  that  the  effects  of  injection 
angle,  coolant  mass  flow  rate,  Mach  number,  upstream  boundary 
layer  thickness,  slot  width,  and  the  presence  of  upstream  cooling 
slots  can  be  investigated. 


McFarland , E . and  Tabakoff , VI.  , "Study  of  Impingement  Heat- 
Transfer  With  Rough  Surfaces,"  Report  No.  75-46,  August  1975 . 

Abstract 

The  effects  of  surface  roughness  on  impingement  cooling  of 
a gas  turbine  blade  were  determined  experimentally.  The  study 
was  conducted  in  a heated  flow  two-dimensional  cascade  tunnel. 
Results  showed  that  roughening  of  the  impingement  surfaces 
decreased  the  cooxing  effectiveness  from  that  of  a smooth 
surface  in  an  optimal  cooling  configuration.  However,  the 
roughened  surfaces  increased  the  cooling  effectiveness  over 
that  of  a smooth  surface  in  a nonoptimal  configuration. 
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FILM  COOLING  EFFECTIVENESS  FOR  COMBUSTION  CHAMBERS 


W.  TABAKOFF  and  R.  RAVURI 
University  of  Cincinnati 
Cincinnati,  Ohio 


Abstract 

The  most  effective  technique  for  cool- 
ing combustion  chamber  walls  is  accomplish- 
ed by  the  utilization  of  a continuous  slot 
jet,  which  produces  a thin  sheet  of  cool 
air  on  the  wall  as  a protection  against 
convective  heating  by  the  hot  gases.  A.i 
alternative  method  of  cooling  is  the  injec- 
tion of  cold  gas  through  discrete  louvers 
cut  into  the  wall.  When  the  louvers  are  cf 
small  width,  the  problem  becomes  similar  to 
film  cooling  by  a row  of  discrete  jets.  In 
this  investigation  the  cooling  effective- 
ness of  a row  of  discrete  jets  is  compared 
experimentally  with  that  of  a continuous 
slot  jet  having  an  equivalent  qoolant  flow 
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Film  Cooling  by  Oblique  Slot  Injection 
in  High-Speed  Laminar  Flow 

R.  H.  Nilson*  and  Y.  G.  Tsueit 
University  of  Cincinnati,  Cincinnati,  Ohio 


Wall-cooling  effectiveness  is  investigated  for  oblique  injection  of  coolant  through  single  or  multiple  wai)  site’s 
into  a high-speed  laminar  compressible  boundary  laser  bv  numerical  solutions  of  the  boundary-laser  etpia'inns. 
A giid  control  procedure  which  maintains  a consiun!  flow  rate  between  grid  lines  is  found  lo  be  well  soiled  to  the 
present  injection  calculations  wherein  the  boundary  -layer  growth  in  the  slot  is  as  much  as  a hundred-fold  and  the 
Ingiludinal  component  of  the  injection  velocity  is  in  some  cases  as  huge  as  the  freestream  velocity  . J ilm-cooling 
effectiveness  is  re  ported  for  a variety  of  injection  configurations'  so  that  the  effects  of  coolant  mass  floss  rate,  in  - 
jedion  angle,  ttpst; earn  boundary -lay  er  thiekn-  ,s,  sbd  width,  and  tire  presence  of  upstream  cooling  slots  can  he 
investigated,  i or  the  coolant  mass  flow  rates, considered,  normal  injection  provides  better  cooling  than  tangen- 
tial injection,  particularly  when  frictional  heating  effects  caused  by  tangential  injection  become  a dominant  con- 
sideration. However,  the  excessive  boundary-layer  growth  which  accompanies  normal  injection  may  reduce 
aerodynamic  performance,  thus  making  inclined  injection  a desirable  compromise. 


Nomenclature 

C_  = constant  pressure  specific  heat 
/.  - length  of  leading  edge  upstream  of  injection  slot 

A/ o,  - freestream  Mach  number 

mc  = coolant  mass  flow  tale 

Pr  = P:  tndtl  number  (taken  as  1 .0  in  present  study) 
s - siot  width 

T - temperature 

Ta  = adiabatic  wall  temperature  upstream  of  injection 

slot 

7\  = coolant  temperature 

Tw  = wall  temperature 

- freestream  temperature 

U , — longitudinal  component  of  injection  velocity 

Ua  = freestream  velocity 

n - longitudinal  velocity  component 

VJ  - transverse  component  of  injection  velocity 

x ,y  = coordinates  parallel  and  normal  to  main  flow 

re  - injection  angle  as  show  n in  Fig.  3. 

9*  = - displacement  thickness  defined  in  Eq.  (6) 

17  = cooling  effectiveness  defined  in  Eq.  (5) 

v = stream  function  defined  in  Eq.  (2) 

v - kinematic  viscosity 

p = density 

r,  = shearing  stress  at  wall 

w = dimensionless  stream  function  defined  in  Eq.  (4) 

Introduction 

WALL  cooling  or  reduction  in  heat  transfer  by 
boundary-layer  injection  is  u-ed  in  engineering  ap- 
plication..1 Although  injection  geometry  varie-,  all  con- 
figurations arc  i;  eluded  in  either  of  tw  o 1 •••gores:  I)  parallel 
injection  in  which  a layer  of  coolant  enters  beneath  the  boun- 
dary layer  through  an  offset  wall;  and  2)  transverse  injection 
where  coolant  is  blown  up  ini  the  boundary  layer  through  a 
slot  o;  holes  in  the  wall.  In  either  of  thc  >c  cases  the  angle  be- 
tween the  coolant  flow  and  the  primary  flow  may  vary  from 
Received  October  7,  1974;  revision  received  February  10,  1975. 
Research  partially  sponsored  hv  U.S.  Army  Research  Of- 
fice—Durham.  under  Contract  DAHC 04-59-0-0016. 

IrJ'x  categories:  Boundary  Layers  and  Convective  Heat  Tran- 
sfer I aminar;  Jets,  Wakes,  ansi  Vrcid-Invireid  Flow  Interaction. 
‘Instructor,  Department  of  Mechanical  Engineering. 

^Associate  Professor,  Department  of  Meehar  cal  Engineering. 
Member  AIAA. 


tangential  to  normal.  Most  previous  investigations  are  con- 
fined either  to  parallel  injection  which  is  tangent  to  the 
primary  flow  or  to  transverse  injection  which  is  normal  to  the 
primary  flow.  The  present  study  considers  the  film-cooling  ef- 
fectiveness of  transverse  injection  at  various  injection  angles 
through  single  or  multiple  wall  slots  into  a compressible 
boundary  layer.  Results  have  been  reported  for  low-speed 
laminar  flow  (U-  50  m/sec).2  High-speed  laminar  flow  is 
reported  herein,  and  work  is  presently  underway  to  extend  the 
study. 

• Numerical  solutions  of  the  boundary-layer  equations  are 
obtained  by  a finite-difference  method  which  has  been  ex- 
tensively tested.  Film  cooling  effectiveness  is  presented  for  a 
variety  of  injection  configurations  so  that  the  effects  of 
freestream  Mach  number,  coolant  mass  flow,  injection  angle, 
boundary-layer  thickness,  siot  width,  and  the  presence  of  up- 
stream cooling  slots  can  be  investigated. 

Numerical  Method 

A finite-difference  method2  is  used  to  solve  the  following 
form  of  the  boundary-layer  equations 


du 
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die 

dp  < 

Plt  ~ 
dx 

~ dx 

dy 

dx  * < 

dT 

pu  

dx 
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it  dp 
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+ 9y 

Pr 

where  the  stream  function  satisfies  the  requirements 
dv  (■•”  9 

T—  =pu  = — [Po^'o—  \ T—  (PM)  dr]  (2) 

dy  9.v  Jo  dx 


The  solution  method  is  quite  general  and  allows  arbitrary 
specification  of  equation  of  state,  viscosity  ntodci,  prewar e 
distribution,  wall  temperature  or  heat  flux  distribution,  a- 
well  as  arbitrary  disttibution  of  both  transverse  and 
longitudinal  velocity  components  at  the  injection  slot. 

A system  of  finite-difference  equations  is  derived  '■> 
integration  of  the  boundary-layer  equ  it  c 
trol  volume  and  after  linearization  the  difference  •....  • 
are  reduced  to  the  tri-diagonal  forms: 
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A Numerical  Solution  for  the  Heat  Transfer 
Between  an  Axi-Sym  metric  Air  Jet  and 
Heated  Plate 
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R.  RAVURI 

Teledyne  CAE 
Toledo,  Ohio 


W.  TABAKOFF 

Department  of  Aerospace  Engineering. 
University  of  Cincinnati. 
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Impinging  air  jets  aie  being  used  extensively  in  the  industry  due  to  the'f  characterist'c  high 
rate  of  heat  transfer.  A numerical  scheme  is  devised  to  predict  the  heat  transfer  between  an 
axi-symme'r  c am  jet  and  a heated  plate,  neglecting  the  effect  of  curvature,  compressibility, 
and  turbulence.  In  this  analysis  the  momentum  and  continuity  equations  are  transformed  to 
vorticity  and  stream  function  equations,  respectively,  and  soivect  by  iterative  successive 
substitution  techniques.  Tne  flow  field  is  determined  which  subsequently  allows  for  a 
calculation  of  the  temperature  distribution  in  the  fluid.  It  is  possible  to  predict  the  average 
heat  transfer  coefficient  of  a row  of  axi-symmelnc  air  jets  impinging  on  a heated  plate  by 
assuming  that  each  jet  cools  a fraction  of  the  plate  without  any  interference  from  the 
neighboring  jets.  The  heat  transfer  coefficients  are  compared  with  some  experimentally 
determined  values. 
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Copies  will  be  available  until  August  1.  1976. 
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The  effects  of  surface  roughness  on  impingement  cooling  of  a gas  turbine  blade  were 
experimentally  determined.  The  study  was  conducted  in  s heated  (low  two-dimensional 
cascade  tunnel.  Results  showed  roughening  impingement  surfaces  decreased  trie  cooling 
effectiveness  from  that  of  a smooth  surface  in  an  optimal  cooling  configuration  (H  /D  = 1 ). 
However,  the  roughened  surfaces  increased  the  cooling  effectiveness  over  that  ot  a smooth 
surface  which  was  not  in  an  optimal  configuration  (H/D  - 2). 
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